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FOREWORD 
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F.  J.  Smith,  LMSC  Program  Engineers.  The  suggestions  by  Dr.  J.  F.  Riley  are 
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This  technical  report  has  been  reviewed  and  approved. 


„  W-'  ^  _ 

FRANKLIN  C.  DAVIES,  Colonel,  USAF 
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ABSTRACT 


The  objective  of  this  research  was  to  study  the  incendiary  potential  of  inter- 
metallic  reactions.  A  thorough  literature  review  of  the  thermochemistry  of  exothermic 
alloy  formation,  measurement  of  the  reaction  rates,  reaction  temperatures,  and 
ignition  characteristics  of  loose  and  pressed  binary  and  tertiary  mixed  powders 
provided  basic  information  for  incendiary  assessment.  The  borides,  carbides,  and 
aluminides  of  titanium,  zirconium,  and  nickel  were  found  to  offer  the  greatest  promise. 
The  data  were  compared  with  those  obtained  by  a  computerized  heat  transfer  analysis 
of  rod  shaped  and  hollow  cylindrical  geometries.  The  rate  controlling  parameters 
were  found  to  be  particle  size,  thermal  conductivity,  heat  of  reaction,  and  diffusivity. 
In  most  cases  it  was  found  that  reaction  characteristics  were  independent  of  particle 
size  distribution,  particle  shapes,  contamination,  bulk  density,  and  compaction. 
Experimental  data  agreed  with  analytical  predictions.  Interm etallic  reactions  were 
deemed  to  be  a  novel  class  of  reactions  and  the  results  were  found  to  have  wide  appli¬ 
cability  in  improving  the  performance  characteristics  of  pyrotechnic,  ordnance,  and 
incendiary  devices. 


Distribution  limited  to  U.S.  Government  agencies  only;  this  report  documents "th'e 
study  of  the  incendiary  potential  of  intermetallic  reactions;  distribution  limitation 
applied  July  1971.  Other  requests  for  this  document  must  be  referred  to  the  Air 
Force  Armament  Laboratory  (DLIF),  Eglin  Air  Force  Base,  Florida  32542. 
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SECTION  I 


INTRODUCTION 


Pyrotechnic  reactions  are  of  two  kinds:  (1)  gas  phase  reactions  which  are 
controlled  by  pressure  and  surface  area,  such  as  black  powder  and  single  or  double 
base  propellants;  and  (2)  condensed  phase  reactions  which  are  controlled  by  thermal 
and  molecular  transport,  such  as  thermite  and  intermetallic  reactions. 

Exothermic  intermetallic  reactions  are  of  practical  importance  because  of  their 
high  heat  of  reaction  and  the  magnitude  of  the  resulting  temperature.  Intermetallic 
reactions  occur  with  the  release  of  a  heat  of  solution  as  the  result  of  the  reaction  of 
finely  powdered  metals.  If  the  heat  released  in  reacting  or  forming  an  alloy  is  suf¬ 
ficient,  the  reaction  will  provide  an  intense  source  of  thermal  radiation.  The  result¬ 
ing  temperature  may  be  sufficient  to  melt  or  even  vaporize  some  constituents  which  can 
then  react  with  oxygen  .releasing  additional  heat.  Intermetallic  reactions  which  provide 
intense  sources  of  thermal  radiation  and  high  temperatures  have  potential  applications 
in  advanced  incendiaries  and  pyrotechnics. 

One  purpose  of  this  study  was  to  review  the  thermochemistry  of  exothermic 
intermetallic  reactions  and  to  present  a  discussion  of  such  condensed  phase  processes. 
In  keeping  with  the  generally  accepted  scope  of  intermetallic  reactions,  boron,  carbon, 
and  silicon  have  been  considered  to  be  metallic  while  phosphorus,  sulfur,  arsenic,  and 
selenium  have  been  excluded. 

Other  objectives  of  this  research  were  to  measure  the  heats  of  reaction,  the 
reaction  temperatures,  reaction  rates.and  ignition  characteristics  for  promising 
candidates;  to  determine  the  effects  of  composition,  particle  size,  compaction,  and 
reaction  mechanism  on  the  incendiary  potential;  to  compare  these  data  with  those 
obtained  on  thermate;  and  to  provide  a  scientific  and  engineering  basis  for  expla'  T  ng 
the  findings. 

The  term  "incendiary  potential"  refers  to  thermochemical  and  kinetic  character¬ 
istics  and  docs  not  include  the  assessment  of  the  effects  of  these  reactions  on  a  target. 
As  a  result  of  this  study,  intermetallic  reactions  were  found  to  provide  a  novel  means 
of  improving  the  reliability  and  reaction  characteristics  of  existing  thermate-based 
incendiaries  and,  furthermore,  to  alter  ignition  delays  and  ignition  energy  requirements 
of  other  types  of  pyrotechnic  devices  such  as  smoke  cartridges,  signal  markers,  delay 
fuzes,  and  tracer  munitions. 
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SECTION  n 

THE  THERMOCHEMISTRY  OF  INTERM ETALLIC  REACTIONS 


1.  LITERATURE  SURVEY 

A  thorough  review  of  the  thermochemistry  of  exothermic  intermetallic  reactions 
was  undertaken  to  assure  complete  coverage  of  all  promising  binary  combinations  of 
elementary  constituents.  The  search  focused  on  two  areas:  equilibrium  thermodynamic 
data  on  the  heats  of  fonnation  (heats  of  reaction),  and  applications  of  exothermic  inter¬ 
metallic  reactions  to  problems  of  a  metallurgical  or  pyrotechnic  nature. 

Thermodynamic  data  are  plentiful.  Unfortunately,  little  agreement  exists  and 
rarely  do  the  authors  define  their  reference  conditions  with  care.  The  laborious 
task  of  critical  evaluation  of  the  published  data  has  been  performed  during  the  last 
two  decades  by  R.  Hultgren  and  his  associates  at  the  University  of  California  (1,2). 

As  a  result  of  the  study  presented  here,  it  was  found  that  favorable  thermo¬ 
dynamics  is  a  necessary  but  not  a  sufficient  condition  to  insure  a  self-sustained 
reaction.  Special  experimental  and  analytical  studies  were  performed  to  define  the 
additionally  required  conditions. 

In  the  survey,  a  few  applications  of  exothermic  intermetallic  reactions  were 
found  such  as  a  commercial  Al-Pd  composite  (3)  called  "Pyrofuze"  and  a  reported 
similar  application  of  magnesium -platinum  and  lead-lithium  (4).  More  recent 
applications  are  discussed  in  Section  V  of  this  report. 

Tables  I  through  VIII  list  the  exothermic  heats  of  reaction  of  binary  mixtures 
of  various  elements  of  the  first  through  the  fourth  groups  of  the  periodic  table  with 
elements  of  the  third  through  the  eighth  groups.  The  compounds  listed  under  the 
heading  "System"  indicate  the  stoichiometric  proportion  of  the  reagents  as  well  as 
the  end  product  to  which  the  thermal  data  have  reference.  It  must  be  understood 
that  the  reacting  system  consists  of  mixtures  of  the  finely  powdered  elements,  and 
the  conclusions  drawn  in  this  study  have  validity  regardless  of  whether  the  actual 
phase  described  by  the  stoichiometric  formula  is  formed. 

Heats  of  formation  are  reported  on  a  molar  basis  for  the  compound  as  written. 
They  are  equal  to  the  heats  of  reaction  where  the  starting  materials  are  elements 
in  their  standard  state.  The  literature  source  cited  refers  to  this  enthalpy  value. 
Assuming  that  a  system  consists  of  a  stoichiometric  blend  of  the  powders  compressed 
to  the  bulk  density  of  the  parent  elements,  it  is  possible  to  compute  the  energy  density 
as  well.  The  density  of  the  intermetallic  product  is  usually  significantly  greater 
than  the  density  of  the  unreacted  mixture.  Because  the  theoretical  bulk  density  is 
rarely  achieved,  intermetallic  reactions  can  be  expected  to  generate  a  product  with 
an  appreciable  void  fraction. 
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TABLE  I.  INTERMETALLIC  REACTIONS  OF  THE  ALKALI  METALS 


TT 


J0 


*4 


f 

■'  K 


.*• 


tr. 


¥ 

i- 


£ 

fr 


.js.^uai 


*■ 


Heat  of  Reaction 

Adiabatic 

Reaction 

Temperature 

<°C>* 

Lowest 

Melting 

Point 

<°C) 

Group 

System 

kcal/mole 

cal/ g 

cal/cm3 

Ref. 

m 

LiAl 

-15.1 

-443 

-555 

860  (1) 

109 

nn 

rv 

LiSn 

-16.8 

-134 

-575 

1000  (1) 

109 

sn 

LiPb 

-14.6 

-68 

-465 

960  (1) 

109 

■pH 

NaSn 

-10.2 

-71.4 

-254 

800  (s) 

98 

HI 

Na2C2 

-  9.6 

-137 

-161 

460  (s) 

98 

1111 

V 

LiSb 

-21.8 

-170 

-700 

1060  (1) 

109 

WM 

LiBi 

-18.4 

-85.5 

-535 

1000  (1) 

109 

NaSb 

-15.8 

-109 

-375 

780  (1) 

98 

i 

KBi 

-18.0 

-55 

-204 

9R0  (1) 

64 

i 

*  Based  on  an  ambient  temperature  of  298 °K. 


The  adiabatic  reaction  temperature  T„  is  that  temperature  at  which  the 
enthalpy  of  formation,  A  Hf ,  equals  the  enthalpy  of  the  product,  namely 

(A  Hf)T^  =  (Ht  -  HgT)T^ 

Where  the  enthalpies  of  the  intermetallic  product  are  not  known,  the  enthalpies 
of  the  starting  materials  were  used  (1,2): 

(A  Hf)T^  =  (Ht  -  «st)a>Tr  +  (ht  "  Hs^B,  Tr 

The  state  of  the  product  at  the  adiabatic  reaction  temperature  was  checked 
against  the  phase  diagram  and  the  appropriate  symbol  (s,  s-1,  1,  and  1-v)  was  used 
to  denote  solid,  liquid,  and  vapor  states.  The  adiabatic  temperature  is  the  maximum 
temperature  which  the  reaction  can  achieve.  It  is  the  upper  boundary  condition 
imposed  upon  a  computer  code  calculation  of  intermetallic  reactions  as  well  as  a 
guide  in  predicting  whether  a  self-sustained  intern.  ,;tallic  reaction  is  possible. 

The  initial  step  in  the  reaction  process  is  thought  to  be  solid-solid  diffusion 
between  the  reagent  particles.  The  resulting  heat  of  reaction  will  not  be  conducted 
away  from  the  reaction  zone  if  the  thermo!  conductivity  of  the  mixture  is  low. 
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TABLE  II.  INTERMETALLIC  REACTIONS  OF  THE  ALKALINE  EARTHS 


Group 

System 

Heats  of  Reaction 

Adiabatic 

Reaction 

Temperature 

(°C)(*) 

Lowest 

Melting 

Point 

<*C) 

Ref. 

kcal/mole 

H 

cal/cm^ 

I 

LiMg 

-2.56 

109 

1 

CuMg 

-5.0 

650 

2 

II 

CaMg9 

-9.4 

-106 

-175 

540  (s) 

G50 

2 

III 

Al9Ca 

-52.5 

-560 

-1150 

1440  (1-v) 

650 

i 

YMg 

-31.0 

-274 

-912 

1670  (1) 

650 

7 

CeMg 

-18.3 

-111 

-523 

920  (1) 

650 

7 

U  Mg 

-41.4 

-157 

-1570 

1940  (1) 

650 

7 

UBe13 

-39.3 

-111 

-520 

500  (s) 

1132 

8,2 

PuBc j3 

-35.7 

-100 

-476 

450  (s) 

640 

2 

ThMg2 

-15.2 

-54 

-317 

700  (s) 

650 

9 

IV 

SnMg., 

-18.9 

-113 

-450 

890  (s) 

890 

1 

SnCa9 

-75.0 

-379 

-1110 

1484  (1-v) 

232 

1 

SnBa9 

-90.0 

-240 

-1000 

2125  (1-v) 

232 

1 

PbCa9 

-50. 7 

-172 

-705 

1440  (1-v) 

327 

1 

PbBa9 

-69 

-143 

-715 

1640  (1-v) 

327 

1 

v 

NbBe,- 

-46.4 

-336 

-1300 

1390  (s) 

1277 

10 

Sb9Ba9 

-87.5 

-133 

-563 

1560  (1) 

630 

1 

Bi9Ba 
_ I _ _ 

-72.5 

-87.3 

-169 

1400  (1) 

271 

1 

(*)  Based  on  an  ambient  temperature  of  298 °K. 


Because  of  the  small  size  of  the  individual  grains,  the  resulting  temperature  rise 
may  therefore  be  sufficient  to  cause  the  reactant  or  the  product  to  melt  so  as  to 
destroy  the  diffusion  barrier  and  expose  fresh  reagent.  If  no  melting  occurs,  the 
diffusion  barrier  may  reduce  the  reaction  rate  to  the  point  of  quenching.  The  most 
important  factors  which  promote  quenching  are; 
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INTERMETALLIC  REACTIONS  OF  GROUP  IU  ELEMENTS:  BORIDES 


Heats  of  Reaction 

Adiabatic 

Lowest 

Reaction 

Temperature 

Melting 

Point 

Group 

System 

kcal/mole 

cal/g 

cal/cm^ 

Ref. 

- 

. 

(*C)(*) 

(°C) 

n 

MgB6 

-22.4 

-251 

-600 

645  (s) 

650 

11 

r  - 

in 

YB6 

-24 

-156 

-503 

700  (s) 

1509 

12,13 

•  1 

LaBg 

-114 

560 

-2410 

2230  (1-s) 

1509 

12,13 

■ 

IV 

CB4 

-17 

-326 

-758 

750  (s) 

2030 

14 

’)* 

SiB„ 

-7 

-76.4 

-177 

230  (s) 

1410 

15 

Jl 

-71.6 

-1155 

-4040 

3770  (1) 

1668 

2 

ZrB2 

-78. 0 

-690 

-3330 

3400  (1) 

1852 

2 

: 

HfB2 

-78.6 

-394 

-3450 

3380  (s-1) 

2030 

2 

• 

V 

VB 

-32.6 

-528 

-2500 

2400  (s) 

1900 

16 

NbB2 

-58.6 

-511 

-2870 

2900  (s) 

2030 

17 

TaB^ 

-43.9 

-216 

-2610 

2400  (s) 

2030 

17 

$ 

VI 

CrB2 

-30 

-408 

-1880 

1500  (s) 

1875 

13 

MoB2 

-23 

-196 

-1280 

1260  (s) 

2030 

13 

W0B 

—  35 

—  83 

-1350 

960  (S) 

2030 

13 

r»*> 

*  ■ 

-  D 

1 

VII 

MnB2 

-19 

-248 

-1140 

800  (s) 

1245 

13 

(*)Based  on  an  ambient  temperature  of  298°K. 


(1)  High  thermal  conductivity 

( 2)  Large  particle  size 

(3)  Adiabatic  reaction  temperature  too  low  to  promote  adequate  diffusion 

He:  ce,  it  is  concluded  that  where  the  adiabatic  reaction  temperature  does  not  exceed 
the  melting  point  of  the  product,  the  reactions  are  not  self-sustaining.  This  con¬ 
dition  is  shown  by  the  notations)  behind  the  adiabatic  reaction  temperature.  Examples 
are  CaMgg,  UBei3,  ThMg2,  NbBes,  MgB6>  YB6,  CrB2.  M0B2,  W2B5,  MnB2,  CuAl, 
CrAl3,  MnAl,  SmB^,  CaC2.  B4C,  VC,  NbGe2>  NbC,  TaC,  Cr?C3,  Mo2C ,  WC, 
Mn7C3>  Mg2Si,  Cr5Si3,  MnSi2,  FeSi,  CoSi,  NiSi,  WSi2,  and  CuPd.  Many  of  the 
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TABLE  IV.  INTERM ETALLIC  REACTIONS  OF  GROUP  HI  ELEMENTS:  ALUMINIDES 


0 

jft 

rl* 

% 


■  # 

JL'. 


0 


II 


f 


ff 


'  m 


h 


Group 

System 

- -  — T 

Heats  of  Reaction  ! 

Adiabatic 

Reaction 

Temperature 

(°C)<*> 

Lowest 

Melting 

Point 

(°C) 

Ref. 

kcal/mole 

cal/g 

cal/cm^ 

I 

CuAl 

-9.8 

-756 

870  (s) 

18 

LiAl 

-15.1 

-443 

-555 

860  (1) 

5 

II 

CaAl2 

-52.5 

-5G0 

-1 150 

1440  (1-v) 

1,18 

III 

B2A1 

-3G 

-1850 

1900  (s) 

G60 

19 

LaAl . 

4 

-41.2 

-1G6 

-780 

1222  (s-1) 

GG0 

1 

IV 

TiAl.^ 

-28.8 

-314 

-1100 

1370  (1-s) 

6G0 

2 

ZrAl2 

-40.8 

-281 

-1200 

1650  (1-s) 

660 

20 

V 

vai3 

-20. 0 

-198 

-792 

750  (s) 

660 

18 

VI 

CrAl3 

-JG 

-120 

-430 

520  (s) 

660 

1,18 

vn 

MnAl 

-10.2 

-124 

-58G 

530  (s) 

G60 

1,18 

VIII 

FeAl 

-12.14 

-147 

-725 

1150  (s) 

660 

1,19 

CoAl 

-25.4 

-308 

-IG00 

158C  (s) 

660 

1 

NiAl 

-28.1 

-329 

-1710 

1G39  (s-1) 

660 

2 

PdAl 

-43.7 

-327 

-2890 

2380  (1) 

660 

3 

PtAl 

■  -48 

-21G 

2510 

2800  (1) 

660 

2 

(*)Based  on  an  Ambient  Temperature  of  298°K. 


other  systems  are  marginal  in  that  iow  heat  output  would  make  a  purely  incendiary 
application  questionable  or  high  thermal  conductivities  would  cause  quenching. 3 
Systems  with  the  maximum  potential  (energy  density  in  excess  of  2000  cai/cm 
and  reaction  temperature  in  excess  of  2000°C)  include  LaB(j,  TiB2 ,  ZrB2>  PdAl, 
PtAl,  BeC2>  TiC,  ZrC,  VSig,  and  NbSia.  Because  of  high  cost  and  short  supply, 
one  should  exclude  systems  requiring  La,  Pd,  and  Pt.  The  systems  selected  on 
this  basis  permit  a  study  over  a  wide  range  of  compositions,  particle  size,  and 
ambient  temperature  and  resulted  in  sevoial  new  concepts  for  pyrotechnic  or 
incendiary  applications. 
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TABLE  V.  1NT ERMETALLIC  REACTIONS  OF  GROUP  ILL  ELEMENTS:  HAKE  EARTHS 
AND  ACTINIDES 


Group 

System 

Heats  of  Reaction 

Adiabatic 

Reaction 

Temperature 

(°C)(*) 

Lowe  st 
Melting 
Point 
(°C) 

Ref. 

kcal/mole 

cal/g 

cal/cm3 

Rare 

SmB6 

-50 

-232 

-1050 

960  (s) 

1072 

22 

Earths 

CeB6 

-84 

-410 

-1850 

1700  (s) 

795 

12 

CeAl4 

-31.2 

-126 

-458 

900  (s) 

660 

1 

PrAl4 

-54.0 

-216 

-800 

1430  (s-1) 

660 

1 

Acti- 

ThB 

-52.0 

-189 

-1360 

1550 (s) 

1750 

12 

nides 

ub4 

-60 

-214 

-1950 

1770  (s) 

1132 

12 

uai4 

-29.8 

-86 

-565 

760  (s) 

660 

23 

PuA14 

-43.2 

-123 

-820 

1130  (s) 

G60 

24 

l(*) Based  on  an  ambient  temperature  of  298  °K. 


2.  EXPERIMENTAL  STUDIES 
a.  Materials  and  Equipment 

Starting  materials  were  procured  from  commercial  sources  and  used  without 
any  further  attempt  at  purification  or  sieving.  Appendix  I  lists  the  sup¬ 
pliers  and  the  nominal  specification  of  these  powders.  The  powders  were  further 
characterized  by  scanning  electron  microscopy  which  permitted  a  determination 
of  geometry,  size,  and  size  distribution  through  measurements  of  the  resulting 
photographs.  Appendix  I  also  lists  the  most  typical  or  rate  determining  particle 
size.  It  is  this  value  which  will  be  used  in  the  subsequent  discussion  of  the  data. 

Some  powders  are  sensitive  to  air  and  are  shipped  under  oil  or  water;  these 
were  dried  in  a  vacuum  chamber  and  manipulated  in  an  argon-filled  glove  box. 

Determinations  of  the  reaction  temperatures  and  emissivities  were  performed 
using  fast  response  pyromotry.  For  events  of  brief  duration,  the  calibrated  anode 
current  of  a  photomultiplier  tube  was  used  to  measure  brightness-temperature.  The 
response  time  of  this  signal  was  a  fraction  of  a  millisecond,  and  the  signal  was 
monitored  using  a  fast  oscilloscope  or  a  recording  oscillograph/  Two  instruments 


TABLE  VI.  INTERMETALLIC  REACTIONS  OF  GROUP  IV  ELEMENTS:  CARBIDES 


Group 

System 

Heats  of  Reaction 

Adiabatic 

Reaction 

Temperature 

(“C)H 

Lowest 

Melting 

Point 

(°C) 

Ref. 

kcal/mole 

cal/g 

cal/  cm^ 

I 

Na2C2 

-9. 6G0 

138 

196 

410  (s) 

98 

6 

II 

BeC2 

-57.7 

-1750 

-3720 

2770  (1-v) 

1277 

13 

CaC2 

-15.0 

-222 

-388 

900  (s) 

838 

6 

SrC9 

-20.2 

-182 

-460 

1300  (s) 

768 

25 

BaC2 

-19.5 

-121 

-391 

1200  (s) 

714 

25 

ffl 

b4c 

-17 

-326 

-758 

750  (s) 

2030 

14 

ai4c3 

-53.4 

-371 

-965 

1400  (s) 

660 

26 

l.c2 

-29 

-178 

-870 

1700  (3) 

920 

27 

CeC9 

-27.9 

-170 

-880 

1560  (s) 

795 

28 

uc/ 

-25 

-95.5 

-1082 

1300  (s) 

1132 

29 

ThC2 

-46 

-179 

-1280 

2800  (S) 

1750 

29 

IV 

SiC 

-15.8 

-394 

-910 

1680  (s) 

1410 

30 

TiC 

-44.1 

-737 

-3270 

3600  (1) 

1668 

2 

ZrC 

-47.0 

-461 

-2460 

3800  (1) 

1852 

2 

IlfC 

-52.3 

-268 

-2730 

4200  (s) 

2222 

2 

V 

VC 

-25.9 

-412 

-1870 

2030  (s) 

1900 

2 

NbC 

-33.6 

-321 

2730  (s) 

2415 

2 

TaC 

-35 . 5 

-184 

-2380 

2800  (s) 

2896 

2 

VI 

Cr7C3 

-54.50 

-136 

-818 

780  (s) 

1875 

2 

Mo2C 

-11.01 

-57 

-482 

GOO  (s) 

2 

w  c 

-9 . 06 

-49 

-645 

2 

VII 

Mn?C3 

-20.2 

-48 

-300 

300  (s) 

2 

(*)Based  on  an  ambient  temperature  of  298 °K. 


8 


iht 


'  * . 

. 

TABLE  VTL  INTERMETALLIC  REACTIONS  OF  GROUP  IV  ELEMENTS:  SILICIDES 


ft 

Group 

System 

Heats  of  Reaction 

Adiabatic 

Reaction 

Temperature 

(°C)(*) 

Lowest 

Melting 

Point 

‘C 

Ref. 

kcal/ mole 

cal/  g 

cal/cm3 

•:r 

n 

Mg2Si 

-10.2 

-180 

-346 

540  (s) 

650 

31 

i.' 

CaSi 

-36 

-528 

-950 

1670  (1) 

838 

13 

.  f 

in 

Y  Si 

-32.2 

-275 

-1000 

1835  (8-1) 

1410 

13 

CeSi2 

-50.0 

-255 

-1100 

1810  (1) 

795 

13 

ThSi 

-42.0 

-146 

-963 

1590  (s) 

1410 

13 

usi2 

-31.0 

-106 

-940 

1390  (S) 

1132 

32 

•f 

IV 

CSi 

-15.8 

-394 

-910 

1680  (s) 

1410 

30 

•* 

TiSi2 

-32.1 

-308 

-924 

1610  (s) 

1410 

2 

ZrSi2 

-38.1 

-259 

-1020 

1780  (1) 

1410 

2 

V 

vsi2 

-75 

-700 

-2350 

3280  (v-1) 

1410 

13 

NbSi2 

-71.6 

-480 

-2040 

3300  (1) 

1410 

13 

A'  . 

TaSi2 

-27.8 

-117 

-806 

1570  (s) 

1410 

2 

VI 

Cr5Si3 

-77.6 

-227 

-1110 

1230  (S) 

1410 

13 

? 

3 

MoSi? 

-31.3 

-206 

-940 

1550  (s) 

1410 

2 

~  ? 

■>  > 

WSi2 

-22.4 

-93.4 

-670 

1230  (s) 

1410 

13 

i 

VII 

MnSi1  7 

-23.2 

-226 

-847 

1160  (s-1) 

1245 

33 

V1H 

FeSi 

-19.2 

-229 

-1000 

1410  (s-1) 

1410 

13 

CoSi 

-24.0 

-276 

-1285 

1460  (s-1) 

1410 

13 

NiSi 

-20.5 

-236 

-110 

990  (s-1) 

1410 

13,19 

(*)  Based  on  an  am  temperature  of  298°K. 


«»*< 


f 

!»■ 

* 


were  used;  The  Br 
EP-1,  and  the  Infr; 
against  a  Bureau  oi 
against  the  melting  j 


gi  nee  ring  Company's  Photo-Electronic  Pyrometer,  Model 
istrios  Thermodot  Model  TD-llS.  Calibration  was 
ds  lamp.  Temperatures  above  2200°C  were  checked 
f  refractory  metals. 
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TABLE  VHI.  INTERMETALLIC  REACTIONS:  MISCELLANEOUS  SYSTEMS 


Group 

System 

Heats  of  Reaction 

- - 

Adiabatic 

Reaction 

Temperature 

(°C)(*) 

Lowest 

Melting 

Point 

(‘C) 

Ref. 

kcal/mole 

cal/g 

cal/cm3 

I 

NaSb 

-16.2 

-112 

-384 

892  (1-v) 

98 

34 

CuPd 

-7.57 

-44.5 

-472 

600  (s) 

10S3 

2 

n 

My^Gc 

-25.1 

-207 

-605 

1115  (s-1) 

650 

35 

Mg3Sb2 

-38.25 

-121 

-484 

1160  (s) 

630 

2 

MggSn 

-10. 1 

-62 

-145 

500  (s) 

232 

36 

m 

CeZn 

-21.2 

-78 

-535 

1140  (1) 

420 

37,  38 

PuCd 1 , 

-45.9 

-31 

-293 

570  (s) 

321 

37 

PuZn 

-73 . 3 

-71.3 

-600 

700  (s) 

420 

37 

U2Zn17 

-95 . 4 

-60 

-530 

700  (s) 

420 

39 

IV 

PdSn 

-25.24 

-112 

-1000 

1326  (s-1) 

232 

40 

Ce2Pb 

-39.64 

-81.5 

-G65 

1380  (s-1) 

327 

2 

NbGe2 

-20.8 

-84.4 

-524 

1170  (s) 

937 

41 

ZrZn2 

-37.63 

-170 

-1165 

1450  (l,s) 

692 

2 

V 

NbNi 

-10.8 

71 

610 

810  (s) 

1726 

2 

(*)  Based  on  an  ambient  temperature  of  298 °K. 


Assuming  a  groybody  behavior,  the  true  temperature  was  measured  using  a 
Milietron  Two  Color  Pyrometer.  The  response  time  of  this  instrument  was  more 
than  300  msec,  and  it  was  only  of  limited  value  in  the  measurement  of  fast  events. 
Furthermore,  the  assumption  of  a  groybody  was  found  to  be  satisfactory  only  in  the 
study  of  carbide  reactions. 

For  comparative  studies  of  incendiary  effects  and  for  measuring  heat  flux, 
radiometers  were  found  to  be  more  powerful  tools  than  pyrometers.  For  determining 
temperature  profiles  as  a  measure  of  rate  in  fast  reactions  and  a  measure  of  the 
effect  of  mixture  characteristics  on  ignition  characteristics,  fast  response  pyro¬ 
meters  remain  the  preferred  tool. 
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Heats  of  reaction  were  determined  as  a  function  of  composition  in  air  and/ 
or  argon  in  a  Parr  bomb  calorimeter  operating  at  1  atm.  After  blending  and  press¬ 
ing,  the  starting  materials  were  characterized  with  respect  to  uniformity  and  void 
formations  by  sectioning  and  metallographic  examination.  Following  reaction,  the 
degree  of  reaction  was  similarly  characterized. 

Ignition  temperatures  were  measured  by  the  insertion  of  a  fine  thermocouple 
into  the  pressed  reaction  body.  Mixtures  were  initiated  by  wrapping  a  nichrome 
wire  around  the  pressed  specimen  and  heating  the  system  electrically.  It  will  be 
shown  in  this  report  that  ignition  temperatures  are  affected  partly  by  heat  flux  and 
partly  by  mixture  characteristics  such  as  particle  size  and  thermal  conductivity. 

Other  experimental  techniques  relating  to  the  measurement  of  reaction  rates 
and  ignition  rates  will  be  discussed  in  Sections  HI  and  IV  of  this  report. 

Measurements  of  thermal  transport  properties  were  of  importance  in  analyzing 
observed  reaction  phenomena,  and  these  were  made  part  of  this  research  program. 

The  thermal  conductivities  were  determined  from  published  specific  heats  and  measured 
thermal  diffusivities  and  densities.  Thermal  diffusivities  were  measured  using  the 
flash  technique  (Ref.  21)  in  which  a  short  duration  (  <  1  msec)  flash  from  a  xenon  lamp 
is  deposited  on  the  face  of  a  thin  layer  or  wafer  of  intermetallic  mixture,  and  the  time 
required  for  the  back  face  temperature  to  reach  50  percent  of  its  maximum  value  is 
measured  oscillographicaily.  This  value  is  used  in  an  analytical  expression  for  the 
flow  of  heat  through  an  infinite  slab.  Corrections  for  heat  losses  are  applied. 

A  nominal  compaction  pressure  of  20,000  psi  was  used  throughout  this  study 
although  it  was  found  tl  t  intermetallic  reactions  are  not  limited  by  the  compaction 
pressure,  and  that  raos  mixtures  tend  to  expand  to  varying  degrees  when  released 
from  the  die. 

The  performance  of  intermetallic  reactions  was  compared  against  that  of 
thermate,  TH-3,  obtained  for  this  purpose  from  Eglin  Air  Force  Base, 

b.  Experimental  Results 

Table  IX  presents  a  summary  of  the  experimental  data  obtained  in  this  research 
program.  The  measurements  do  not,  in  all  cases,  agree  with  published  data.  This 
is  probably  due  to: 

«  Unreliable  published  data,  particularly  where  these  are  derived  from  solution 
calorimetric  or  electrochemical  data 

•  Incomplete  reaction  because  of  surface  cooling  or  the  presence  of  local 
inhomogeneities 

•  Incomplete  reaction  because  of  dead  sintering 

•  Low  total  enthalpy  changes  lead  to  low  precision  in  calorimetry 

•  Deviations  from  greybody  behavior 
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TABLE  IX.  TABULATION  OF  EXPERIMENTAL  DATA 


System 


Experimental  Heat 
of  Reaction  (cai/j;) 


No  Hx  |  No  Hx 


Measured 

Density 


Hx  No  Rx 


Thermato 


(a)  Computed, 


Measured 

Einissi- 

vlty 

Ignition 
Temp, 
in  "C  Air 

Average 
Reaction 
Temp.  "C 

Reference 
Table  No. 
(This  Report) 

170 

I 

0.  179 

200 

No  publ.  data 

~1 200” 

VI 

0.444 

3010 

0.727 

550 

3000 

in 

530 

2340 

3410 

m 

0.72 

410 

2315 

m 

0.458 

2160 

m 

400 

1200 

IV 

100 

1350 

rv 

640 

1500 

No  publ,  data 

0.376 

640 

1700 

IV 

1100 

600 

2050 

VI 

550 

1540 

VI 

650 

1250 

VI 

450 

1625 

vn 

0. 168 

250 

1430 

vn 

No  publ.  data 

540 

930 

No  publ.  dala 

(47) 

200 

1090 

No  publ.  data 

1820 

I 

200 

1195 

I 

50 

2300 

vrn 

50 

1260 

vm 

50 

1270 

a 

0.275 

600 

2940 

<a) 

The  corresponding  published  data  may  be  found  in  the  earlier  tables  of  this 
section  by  reference  to  the  table  number  given  in  the  right-hand  column  of  Table  IX. 

The  observed  reaction  temperatures  were  in  reasonable  agreement  with  the 
adiabatic  reaction  temperatures.  Wh  re  the  measured  values  were  higher,  it  was 
duo  to  the  ignition  technique  which  involved  heating  the  entire  sample  to  the  ignition 
temperature  so  that  the  reaction  energy  was  augmented  by  the  ignition  energy,  hi 
some  instances,  the  reaction  energy  was  increased  by  reaction  with  air. 

Based  on  the  thermal  measurements  of  this  study,  summarized  in  Table  IX, 
it  was  concluded  that  the  following  systems  have  the  greatest  incendiary  potential: 

•  Boron  with  titanium ,  lithium,  or  zirconium 

•  Carbon  with  titanium  or  zirconium 

•  Aluminum  with  nickel  or  with  mixtures  of  the  above 

This  appraisal,  based  on  the  thermal  data  only,  is  in  general  agreement  with  ' 
the  findings  of  the  literature  survey. 

In  addition  to  these  binary  systems,  the  effect  of  a  third  component  was  experi¬ 
mentally  investigated.  Admixture  of  a  third  component,  such  as  aluminum  to  a  boron- 
titanium  mixture,  results  in  heats  of  reaction  between  those  of  Al-B  and  Ti-B. 

Other  effects  on  the  incendiary  potential  are  due  to  changes  in  thermal  conductivity, 
melting  point  of  an  intermediate  phase,  and  similar  reaction  kinetic  parameters. 

These  effects  are  discussed  in  further  detail  later  in  this  report.  It  was  shown  that 
ease  of  ignition  and  reliability  of  reaction  were  improved  by  lowering  the  thermal 
conductivity.  This  was  demonstrated  by  the  addition  of  small  amounts  of  boron  or  . 
lamp  black  to  reaction  mixtures  of  relatively  high  thermal  conductivity,  such  as 
Ni-Al  or  thermite.  The  low  mechanical  strength  of  some  formulations  may  be 
overcome,  in  part,  by  the  admixture  of  a  binder,  such  as  1  percent  gum  arabic, 
without  noticeable  effect  on  performance  characteristics. 

c.  Analysis  of  Reaction  Products 

Insight  into  the  reaction  mechanism  was  obtained  by  pyromotry  and  by  analysis 
of  the  intermetallie  systems  both  before  and  after  reaction.  Figure  1  shows  a 
section  through  a  pellet  of  Ti-2B  which  is  quite  representative  of  all  the  compressed 
reaction  mixtures.  The  mass  appears  to  be  without  obvious  voids,  the  titanium  (the 
discontinuous  phase)  appears  to  be  uniformly  distributed,  and  the  particles  appear 
to  be  elongated  or  flat  and  preferentially  oriented  (normal  to  the  direction 
of  the  compression)  . 

The  appearance  of  the  reaction  product  reveals  whether  the  reaction  went  to 
completion  and  whether  the  reaction  occurred  in  the  gas  phase  or  in  the  liquid  phase. 
Figure  2  shows  particles  of  titanium  in  a  carbon  matrix  which  are  surrounded  by 
crusts  of  the  reaction  product.  This  layer  of  product  effectively  prevents  further 
reactions.  A  similar  picture  (Figure  3)  shows  particles  of  tantalum  surrounded  by 


Figure  1.  Photomicrograph  of  Typical  Starting  Mixture.  100  x 

Titanium  powder,  30  p  effective  particle  size,  in  contin¬ 
uous  matrix  of  boron  powder,  2  p  particle  size,  pressed 
at  20,000  psi. 
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Figure  2.  Photomicrograph  of  a  Partially  Completed  Inturmetallic  Reaction, 
500  x 

Titanium  powder  particles,  30-p  effective  particle  size  in  a  0.  5-p 
graphite  matrix,  pressed  at  20,000  psi,  heated  to  about  800°C  for 
5  minutes.  Metallographic  section  shows  the  reaction  inhibited  by  a 
crust  of  reaction  product  which  surrounds  the  metal  powder.  Such 
behavior  is  common  when  the  particle  size  or  the  thermal  conduc¬ 
tance  is  too  large.  Such  mixtures  will  sinter  but  not  react  in  a 
self-sustained  manner. 
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Figure  3.  Photomicrograph  of  Incomplete  Intormetallic  Reaction,  5000  x 

Tantalum  powder,  4-p  offective  particle  size  in  2-g  silicon  powder, 
pressed  at  20,000  psi  and  heated  to  abouL  t,00°C  for  5  minutes. 
Metal lographic  section  shows  partial  completion  of  the  x*caction  by 
formation  of  successive  layers  of  the  various  silicide  phases.  This 
type  of  reaction  is  typical  when  the  reaction  temperature  is  insuf¬ 
ficient  to  melt  the  reaction  product.  Such  mixtures  usually  cannot 
be  made  to  react  in  a  self-sustained  manner.  At  times,  preheating 
of  the  reaction  mixture  is  effective  provided  that  the  self-quenching 
shown  in  Figure  2  can  be  avoided. 
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several  phases  of  a  silicide  and  immersed  in  a  matrix  of  silicon.  Such  systems  do 
not  react  in  a  self-sustained  manner.  Usually  it  is  necessary  to  reduce  the  particle 
size  (as  for  titanium  carbide)  or  increase  the  ambient  temperature  so  as  to  raise 
the  reaction  temperature  above  the  melting  point  of  the  product  (as  with  the  tantalum 
silicide)  before  the  reaction  becomes  self-propagating.  Figure  4  shows  the  titanium 
diboride  reaction  product  chat  is  formed  when  the  compressed,  unreacted  mixture 
undergoes  essentially  no  preheating,  prior  to  ignition.  The  reaction  temper¬ 
ature  never  approaches  the  predicted  adiabatic  reaction  temperature  but  is 
sufficient  (2500-3000°C)  to  cause  melting  and  crystallization  of  the  product 
from  the  liquid  phase.  The  appearance  of  voids  and  of  a  fairly  uniform 
reaction  are  typical.  When  a  mixture  is  preheated  above  ambient  temperature, 
the  preheat  temperature  increases  the  adiabatic  reaction  temperature  causing 
the  system  to  react  above  the  boiling  point  of  some  of  the  constituents. 

Such  a  phenomenon  is  accompanied  by  deposition  of  crystalline  product  from 
the  gas  phase  as  shown  in  Figure  5. 
d.  Thermal  Transport  Properties 

Table  X  lists  the  thermal  transport  properties  of  the  reaction  mixtures  selected 
for  this  study.  The  heat  capacities  are  those  of  the  starting  materials  at  ambient 
temperature.  The  densities  are  experimental  values  and  reflect  the  apparent  density 
of  the  compressed  samples.  The  experimental  method  is  designed  to  measure  the 
thermal  diffusivity  which  should,  in  principle,  be  independent  of  density.  Table  X 
shows,  however,  that  thermal  diffusivities  of  compressed  and  uncompressed  samples 
differ  by  as  much  as  300  percent,  indicating  that  conductance  in  fine  powders  is  a 
complex  phenomenon.  It  appears  that,  when  particles  are  finely  divided,  surface 
effects  predominate  and  the  properties  of  the  bulk  material  can  no  longer  be  used. 
There  are  no  good  theoretical  approaches  for  estimating  thermal  conductivities  of  fine 
powders.  For  this  reason,  thermal  transport  data  were  obtained  experimentally  for 
interm etaLlic  mixtures  and  for  the  unmixed,  but  powdered,  reactants.  These  data 
can  be  a  valuable  guide  in  estimating  conductivities  of  mixtures  having  compositions 
which  arc  not  listed  by  serving  as  aids  in  interpolation. 
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Figure  4  .  Into  metallic  Reaction  Product  -  Starting  Mixture  at  Room 
Temperature  500  x 

Photomicrograph  of  TiB2  reaction  product.  The  material 
exhibits  voids  with  incipient  crystal  formation.  This  is  the 
typical  appearance  of  a  product  which  reacted  at  a  temperature 
sufficiently  high  to  cause  melting  but  not  to  produce  a  vapor 
phase.  This  stoichiometric  mixture,  pressed  to  20,000  psi, 
was  ignited  with  a  first  fire  and  within  a  few  milliseconds 
heated  itself  to  about  2800 °C. 
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Figure  5.  Inter-metallic  Reaction  Product,  Starting  Mixture  Preheated,  3000  x 

Photomicrograph  of  a  TiC  reaction  product.  The  material  exhibits 
voids,  crystals  which  appear  formed  from  the  liquid  phase,  and 
large  needle-like  crystals  which  are  believed  to  have  been  grown 
from  the  vapor  phase.  This  is  the  typical  appearance  of  a  product 
which  reacted  at  a  temperature  sufficiently  high  to  cause  partial 
vaporization.  Tills  stoichiometric  mixture  of  titanium  powder  in 
lampblack,  pressed  at  20,000  psi,  was  preheated  to  about  800°C 
when  it  suddenly  heated  itself  to  a  temperature  of  about  3000 ’C 
which  was  retained  for  about  2  sec. 
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E  X.  THERMAL  TRANSPORT  PROPERTIES  OF  iNTERMETALUC  REACTION 
MIXTURES 
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Density 

(g/cnw) 

Thermal  Diffualvity 
(cm2/sec) 

Thermal  Conductivity 
(cai/cm  sec  *K) 

0.975 

7.6  x  10"4 

1.2  10'4 

0.025 

7.19  x  10"4 

1.07  x  lo-4 

2.32 

1.87  x  10'3 

0.7  x  10'3 

0.  571 

9.31  x  10"4 

0  64  x  io~* 

3.28 

3. 85  x  10"2 

1.73  x  10"2 

2.70 

3.81  v  10~3 

1.68  x  lo'3 

2. 10 

2.76  x  id'3 

-3 

1.22  x  10J 

2.95 

1.68  x  lo'3 

6.61  x  lo"4 

3.65 

0.0127 

0.  00352 

5.11 

0.220 

0.156 

4.  50 

0.081 

0.053 

3.95 

0.0134 

0.008 

2.82 

0.0028 

0.0013 

2.44 

0.00345 

0.0014 

2.45 

0.0043 

0.0019 

2.42 

O.UOH5 

0.  0039 

3. 61 

0.00035 

0.00423 

3.76 

0.  009  67 

0.  00447 

1.31 

0. 00090 

0.  00029 

4.35 

0.00374 

0.00107 

1.16 

0.0115 

0.00226 

2.66 

0.451 

0.258 

6.  B0 

0. 052t> 

0.  0316 

2. 48 

0.0117 

0.  0058 

2.75 

0.  0201 

- 

0.  0129 

All  percentages  In  wolght  percent.  For  suppliers  of  metal  powdo^aoo  Appendix 
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1.  EXPERIMENTAL  STUDY 

Two  types  of  reaction  rates  can  bo  defined:  (1)  the  linear  reaction  rate,  cm/sec, 
and  (2)  the  mass  reaction  rate,  g/sec-g.  In  tills  study,  only  the  linear  reaction  rate  will 
be  considered  as  it  is  of  greater  importance  in  pyrotechnic  and  incendiary  applications. 
The  linear  reaction  rate  was  measured  on  bars  pressed  from  powder  mixtures,  as  shown 
in  Figure  G.  The  bar  which  measured  i/4  x  1/4  x  2-1/2  in. ,  was  ignited  at  one  end  with 
a  first  fire. 1  The  rate  of  migration  of  the  reaction  front  was  followed  using  either  a 
high-speed  motion  picture  camera  (Fastax,  1000  f/sec)  or  burn  wires  which  were 
placed  on  the  bar.  Both  methods  gave  equivalent  results.  The  motion  pictures  showed 
that  the  reaction  front  was  typically  straight  and  normal  to  the  length  of  the  bar.  No 
flame  migration  along  the  surface  of  the  bar  was  observed  as  is  the  case  with  solid 
propellant  strands;  hence,  no  inhibition  of  the  sides  was  necessary.  The  reactions 
proceeded  smoothly,  producing  a  solid  rod  which  remained  luminous  for  up  to  a  minute. 
The  reactions,  while  in  theory  not  accompanied  by  a  gas  phase,  usually  generated  a 
cloud  of  particulate  matter  which  was  probably  a  result  of  adsorbed  gases,  the  vaporiz¬ 
ation  of  impurities  (oxides,  magnesium)  from  the  boron , or  the  oxidation  of  the  small 
fraction  of  the  reagents  vaporized  in  the  reaction  zone.  This  smoke  was  sufficient  to 
give  erroneous  emissivity  readings,  particularly  for  the  boron  reactions. 

The  high  reaction  rate  of  intormetallic  reactions ,  as  illustrated  in  Figure  7 , 
was  obtained  from  the  output  signal  of  a  fast  pyrometer  and  recorded  with  an  oscillo¬ 
scope.  The  compressed  mixtui'e  was  ignited  by  an  electrically  heated  niehrome  wire 
which,  by  being  wrapped  around  the  sample,  assured  a  fairly  uniform  temperature 
distribution  throughout  the  sample  at  the  time  of  ignition.  The  temperature  rose  to 
the  peak  within  100msec  during  which  time  the  reaction  was  effectively  completed. 

The  residual  reaction  could  not  have  been  significant,  judging  from  the  steep  temper¬ 
ature  drop  which  followed.  The  reaction  time  can  be  prolonged  by  admixture  of 
coarser  particles  but  at  the  expense  of  peak  temperature. 

2.  EXPERIMENTAL  DATA 

Table  XI  lists  experimental  linear  reaction  rates  fur  various  types  of  intormetallic 
reactions.  Data  on  Pyrofuzo  and  thermate  were  included  for  comparison.  The 
greatest  number  of  parameters  were  investigated  on  Ti-2B;  therefore,  generalized 
conclusions  arc  well  substantiated  for  tiiis  system. 


iIL  is  shown  in  Section  IV  of  this  report  that  no  first  fires  are  necessary  for  the 
initiation  of  intormetallic  reactions.  For  convenience,  however,  the  above  study 
utilized  a  B/BaCrOq  mixture. 
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Mold  for  Burning  Rate  Specimens. 

Sample  of  pressed  reaction  mixture  is 
shown  in  foreground. 
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TABLE  XI.  EXPERIMENTAL  LINEAR  REACTION  RATES 


System 

Effective 

Particle 

Size 

Heat  of 
Reaction 
(cai/g; 

Av.  Heat 
Capacity 
<caVjj-*K) 

Density 

(g/cm3) 

Ambient 

Temperature 

<*K> 

Thermal 

ConthicUvily 

(cal/cm-sec-*K) 

Linear 

Reaction  Rate 
(crn/aoc) 

T1-2B 

1  n  Tl 

2iiE 

1155 

0.306 

2.32 

300 

0.0007 

6.44 

10  n  Ti 

1155 

0.  306 

2.32 

300 

0.0007 

8.33 

2  a  B 

30  m  Ti 

2)i  B 

1155 

0.306 

2.32 

300 

200 

0. 0007 

0. 0007 

6.57 

5.  67 

60  is  Ti 

1155 

0.306 

2.32 

300 

0.0007 

4.38 

2p  B 

60  it  Tl 

10  jjB 

1155 

0.306 

2.32 

300 

0.0007 

1.57 

TI-2B 

10  is  Tl 

925 

0.306 

300 

3.59 

with  20% 
reacted  TiB^ 

2  M  B 

10  fx  TiB2 

an  aboi'u 
with  40% 

TiB2 

as  above 

695 

053 

”  355 

0. 55 

Tl-C 

1  li  Ti 

737 

0.226 

3.38 

300 

0.0173 

T7b3  ■* 

0.  5  »  C 

10  a  Ti 

737 

0.226 

3.38 

300 

0.0173 

1. 49 

0.5  m  C 

Ti-C 

10  ,i  Tl 

590 

0.226 

300 

1. 04 

with  20% 
reacted  TiC 

0.5  »  C 

40  p  TiC 

as  above 
with  40% 

TiC 

as  above 

442 

0.226 

500 

0.30 

Zr-213 

3  ti  Zr 

680 

0.171 

2.95 

300 

0.000561 

0.  89 

2  n  B 

10  fi  Zr 

690 

0.171 

300 

0.  15 

2»*  B 

Zr- C 

1  n  Zr 

0.6  jiC 

461 

0.160 

3.65 

300 

0.00352 

1.86 

.  . 

10  m  Zr 

0.6  M  C 

1.02 

Ti-B 

60  M  Ti 

865 

0.27 

2.76 

300 

0.00108 

3.88 

2  m  B 

T1-3B 

60  n  Ti 

1137 

0.34 

2.10 

300 

0.00122 

3.06 

2  p  B 

Tt-2B 

60  n  Tl 

690 

0.262 

2.45 

300 

0.0043 

2.  59 

2  M  B 

with  30%  Ai 

30  m  Al 

as  above 
with  60%  Al 

ns  above 

033 

0.262 

2.42 

300 

0.0085 

0.62 

Nl-Al 

40  m  N1 

30  M  Al 

329 

0.18 

5.11 

300 

500 

0.156 

no  reaction 
3.50 

as  a bo vo 
with  6%  B 

273 

0.  198 

4.  50 

300 

0.054 

1.73 

with  10%  B 

230 

0.217 

3.95 

300 

0.0082 

0.72 

with  25%  B 

100 

0.273 

2.82 

300 

0.0014 

no  reaction 

Pyrofuze 

NA 

327 

0.  17 

9 

300 

0.3 

20  (48) 

Thormalu 

TH-3 

827 

-0.20 

2.48 

300 

0.0117 

0.276 
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Figure  7.  Temperature-Time  History  of  a  Reaction  Between 
Amorphous  Boron  and  -200  +  325  Mesh  Titanium 


a.  Effect  of  particle  Size 

Linear  reaction  rates,  shown  as  the  average  of  several  runs,  decreased  from 
9.44  to  4.  38  cm/sec  linearly  when  fine  boron  was  used  and  when  the  effective  particle 
size  of  titanium  was  increased  from  1  to  60  p  .  When  the  boron  particles  were  coarse, 
approaching  those  of  titanium ,  the  reaction  rate  was  further  reduced  from  4.  38  to  1. 57 
cm/sec. 

b.  Effect  of  Ambient  Temperature 

The  effects  of  preheating  and  precooling  the  reaction  mixture  are  shown  for 
Ni-Al  and  Ti-2B,  respectively.  The  reasons  for  this  behavior  are  discussed  in  detail 
in  a  subsequent  description  of  analytically  obtained  results.  It  shall  suffice  to  sav 
that  the  ambient  temperature  affects  both  the  ignition  energy  and  the  reaction 
temperature  by  causing  variations  in  the  diffusivity  of  the  reactants. 

c.  Effect  of  Heat  of  Reaction 

Other  systems  have  reaction  rates  lower  than  Ti-2B.  To  resolve  whether  this 
difference  is  due  to  an  inherent  property  of  the  reactants  or  due  to  the  lower  reactivity 
caused  by  a  lower  adiabatic  reaction  temperature,  some  tests  were  made  with  the 
reactants  diluted  with  various  amounts  of  spent  reaction  mixture.  The  effect  of  this 
diluent  was  to  reduce  the  heat  of  reaction  of  the  mixture  without  changing  other  properties 
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such  as  heat  capacity,  diffusion  parameters,  and  thermal  conductivity.  Linear 
reaction  rates  decreased  from  8. 33  to  0.  55  cm/sec  as  the  heat  of  reaction  was 
decreased  from  1155  to  695  cal/g.  A  similar  effect  was  demonstrated  on  Ti-C. 

It  was  postulated,  therefore,  that  a  similar  decrease  of  reaction  rate  would  be  observed 
if  very  large  particles  of  reactant  were  added,  as  these  would  act  as  diluent  in  the  re¬ 
action  zone  by  soaking  up  a  large  amount  of  sensible  heat  before  undergoing  reaction. 

cr 

Similarly,  adding  30  percent  aluminum  to  the  mixture  decreased  the  reaction  rate 
from  4,  38  to  2.  59  cm/sec,  chiefly  because  of  the  adverse  effect  on  heat  of  reaction. 

The  reason  that  the  change  was  not  as  striking  as  in  the  above  case  involving  a  diluent 
is  that  aluminum  increases  the  thermal  conductivity.  This  effect  of  thermal  conduc¬ 
tivity  is  also  visible  by  comparing  the  rates  of  Zr-2B  with  those  of  Ti-2B  containing 
40  percent  TiB2-  Both  systems  have  the  same  heat  of  reaction,  but  Ti-2B  has  a 
higher  reaction  rate,  partially  because  of  its  higher  thermal  conductivity.  The 
effect  of  particle  size,  heat  of  reaction,  and  thermal  transport  will  now  be  considered 
analytically. 

3.  ANALYTICAL  STUDY  OF  REACTION  RATES 
a.  The  Mathematical  Model 

The  absence  of  a  gas  phase  greatly  facilitates  an  analytical  study  of  interm etallic 
reactions.  The  problem  reduces  to  one  of  a  simple  heat  transfer  analysis  uncompli¬ 
cated  by  mass  transfer  considerations  or  movement  of  boundaries.  Consider  a  rod  of 
a  reactive  powder  mixture  which  is  heated  at  one  end  by  an  isothermal  heat  source 
until  the  temperature  has  risen  sufficiently  to  raise  the  diffusion  limited  reaction 
rate  to  a  self-sustaining  value.  The  problem  is  one  of  determining  the  time  to 
achieve  ignition  temperature  (ignition  time)  and  the  rate  of  reaction  zone  migration 
along  the  rod  (linear  reaction  rate). 


To  accomplish  this 
of  effective  size  4' 
diffusivity  given 


end,  assume  that  the  mixture  is  composed  of  particles  A  and  B 
,  thermal  conductivity  k  ,  heat  of  reaction  Q,  and  a  molecular 
by 


D  =  D  exp 
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In  each  slice  sxx  the  rate  of  heat  generation  is  related  to  the  diffusivity  through 
Kick's  Law: 


dQx 

dt 


=  Q 


dC 

dt 


(2) 


where  ip  is  the  thickness  of  the  diffusion  barrier  which  increases  from  0  to  and 
dC  is  the  fraction  of  the  reagent,  reacting  in  time  interval  dt  .  The  heat  loss  from 
each  section  Ax  is  given  by 


dQ2 

dt 


d^T 

K  .  2 

dx 


+ 


ter  27rr  A  x  ,„,4 
tti'2  Ax  ( 


(3) 


Convective  heat  loss  is  ignored  in  this  model.  The  difference  between  the  heat  generation 
and  heat  loss  in  each  section  corresponds  to  heat  accumulation  which  is  reflected  in  a 
rise  of  temperature.  A  computer  program  was  written  which  determines  this  temperature 
as  a  function  of  time  and  number  of  length  increments.  A  listing  of  statements  is  given 
in  Appendix  II.  The  input  variables  are,  in  part,  the  particle  size,  thermal  and 
molecular  diffusivity,  heat  of  reaction,  ambient  temperature,  and  heat  capacity. 

Figure  8  shows  a  sample  of  computed  temperature  profiles  for  a  computer 
experiment  utilizing  data  from  the  titanium -carbon  system.  Several  features  of  the 
experimentally  observed  phenomena  are  reflected  in  these  computed  data,  lending 
confidence  to  the  belief  that  the  actual  reaction  mechanism  is  fairly  accurately  described. 
High  reaction  temperatures  and  the  initiation  of  the  reaction  by  heat  source  at  800°K, 
wore  observed.  This  source  temperature  is  far  below  the  melting  point  of  either  con¬ 
stituent;  thus,  the  assumption  of  a  solid-solid  diffusion  limited  process  appears  reason¬ 
able.  As  intermetallic  reactions  are  equally  well  initiated  in  vacuum  or  in  an  inert 
atmosphere,  the  intermediate  action  of  air  or  of  absorbed  reactants  could  not  play  a 
significant  part.  Initiation  occurred  in  less  than  25  msec,  thereby  again  describing 
an  experimentally  observed  phenomenon.  The  computed  reaction  rate  was  high  compared 
with  the  experimental  value,  probably  because  of  the  fact  that  the  effective  difiusivity 
was  overestimated.  Several  reasons  for  this  are  given  in  the  next  section.  The  computed 
reaction  temperature  approaches  the  adiabatic  value.  Deviation  of  the  computed  reaction 
temperature,  3350°K,  from  the  experimentally  observed  value,  2320°K,  probably 
is  the  result  of  a  deviation  from  air  idealized  state  by  the  experimental  system. 

From  the  described  computer  output,  the  following  data  can  be  obtained: 

•  Ignition  Time.  For  each  source  temperature  and  each  exposure  time,  the 
computer  output  gives  the  time  at  which  the  reaction  roaches  a  self-sustaining 
level. 

•  Ignition  Energy.  The  heat  transferred  during  the  ignition  time  into  the  first 
slice  is  the  ignition  energy.  A  sample  calculation  is  shown  in  Appendix  IV. 

•  Ignition  Delay.  The  time  required  to  raise  the  reaction  from  the  ignition 
temperature  to  the  steady  state  reaction  temperature. 

•  Steady  State  Reaction  Temperature.  The  steady  stale  reaction  temperature 
differs  from  the  adiabatic  reaction  temperature  by  the  fact  that  some  heat, 
however  small,  is  lost  from  the  reaction  zone  to  the  environment. 
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Figure  8.  Computed  Reaction  Characteristics  of  Titanium -Carbon  Reaction  Mixture 


•  Linear  Reaction  Rate.  From  the  difference  in  distance  between  temperature 
profiles  and  the  corresponding  time  increments,  a  velocity  of  reaction  front 
migration  (the  linear  reaction  rate)  can  be  computed. 

•  Reaction  Time.  By  dividing  the  thickness  of  the  reaction  mass  by  the  linear 
reaction  rate,  the  reaction  time  can  be  obtained. 

b.  Results  of  Mathematical  Analysis 

The  computer  program  was  designed  for  parametric  studies  of  reaction 
•haracteristics  and  so  serve  as  a  guide  to  the  most  promising  experimental  studies. 

The  capability  of  the  program  to  portray  actual  system  behavior  depends  on  how  closely 
the  program  simulates  the  reaction  mechanism.  Such  a  test  of  the  computerized 
reaction  model  by  use  of  experimentally  determined  properties  is  shown  in  Table  XII. 

For  this  study,  the  diffusion  parameters,  prefix,  and  activation  energy  were  empirically 
adjustable  parameters.  Also  shown  are  typical  published  data  for  these  parameters  as 
obtained  by  a  layer  analysis  of  borodized and  earborized  metallic  specimens  (43,44). 
Some  published  studies  ignore  the  concentration  dependence  of  diffusivity  and  shifting 
of  phase  boundaries  caused  by  density  changes,  thereby  resulting  in  well  known  lack  of 
reproducibility  in  solid  phase  diffusivity  data.  It  is  felt  that  not  too  much  reliance 
should  be  placed  on  published  diffusivitios,  except  as  guides  in  parametric  analysis, 
because  intormotallic  reactions  proceed  by  formation  of  voids  and  of  products  having 
higher  densities  than  the  reactants.  Systems  of  such  complexity  cannot  be  rigorously 
treated  mathematically 

In  comparing  the  data  given  in  Table  XII,  it  is  found  that  computed  reaction 
rates  are  higher  than  experimentally  determined  rates  except  for  the  Ti-2B  system. 

In  this  system  the  experimental  values  are  higher.  One  reason  for  this,  no  doubt, 
is  the  improper  choice  of  diffusivity  values.  The  discrepancies  may  also  result  from 
changes  of  reaction  mechanism  with  temperature  and  changes  in  conductance  due  to  the 
formation  of  voids.  The  effect  of  the  diffusivity  prefix  on  the  rate  of  a  hypothetical 
interm otallic  reaction  is  given  in  Table  XIII. 

A  parametric  study  of  the  above  hypothetical  system  provided  the  data  in  Table 
XFV  on  the  effect  of  activation  energy  on  reaction  rate. 

The  computer  correctly  predicted  self-quenching  of  reactions  in  systems  hav¬ 
ing  high  thermal  conductivity.  Also,  the  trend  expected  by  the  addition  of  third  com  ¬ 
ponents  or  of  diluents  was  predicted  correctly.  For  instance,  while  third  components 
may  increase  thermal  conductance  and  thereby  increase  the  reaction  rate,  this  effect 
may  be  offset  by  the  lowering  of  the  heat  of  reaction  which  causes  lower  diffusivity 
and  possible  quenching.  A  good  example  is  seen  in  the  addition  of  aluminum  to  Ti-2B. 

This  analytical  approach  is  also  useful  for  obtaining  insight  into  the  complex 
interaction  of  many  parameters.  For  instance,  the  quenching  of  the  Ni-Al  reaction 
was  overcome  when  the  ambient  temperature  was  raised  so  as  to  cause  melting  of 
the  reaction  product  and,  similarly,  lowering  of  the  ambient  temperature  in  Ti-2B 
mixtures  increased  the  ignition  delay  until  so  much  heat  had  soaked  into  the  reaction 
mixture  that  the  apparent  reaction  rate  was  increased.  Similar  heat  transfer  analyses 
of  other  pyrotechnic  reactions  should  be  equally  instructive. 
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Except  where  noted,  the  activation  energy  was  assumed  to  be  20,000  cal/mole. 
Input  Parameters  Not  Listed  are  Given  in  Table  XI. 


TABLE  Xin. EFFECT  OF  DIFFUSIVITY  PREFIX  ON  REACTION  RATE 


Diffusivity  Prefix 
(cm -/see  ) 

Ignition  Time 
(msec) 

Ignition  Delay 
(msec) 

Linear  Reaction  Rate 
(cm/sec) 

5 

1 

5 

0.70 

0.5 

2 

10 

0,46 

0.05 

4 

20 

0.20 

TABLE  XIV,  EFFECT  OF  ACTIVATION  ENERGY  ON  REACTION  RATE 


Activation  Energy 
(eal/mole) 

Linear  Reaction  Rate 
(cm  /sec) 

15,000 

0.97 

20,000 

0.69 

25.000 

0.  50 

30,000 

0.34 

35,000 

0. 12 

Several  refinements  are  available  in  the  program  which  are  designed  to  describe 
actual  physical  characteristics  and  reaction  phenomena  more  precisely. 

(1)  Effect  of  Particle  Shape 

Figure  1  and  scanning  electron  micrographs  show  that  metallic  particles  come 
in  different  shapes  -  some  spherical,  others  like  plalelcltos.  The  e  „.puter  program 
used  in  this  study  assumes  that  the  particles  A  and  B  ,  described  on  page  25,  are 
infinite  plates  of  effective  thickness  'l'  .  Another  program  described  spherical 
particles  A  of  diameter  4  lo  have  been  immersed  in  a  continuous  matrix  ol  B  . 

The  numerical  difference  in  output  data  was  not  sufficiently  large  from  those  obtained 
from  the  parallel  plate  model  to  require  a  correction  of  the  program  presented  in 
this  report. 

(2)  Effect  of  particle  Size  Distribution 

The  program  permits  the  specification  of  15  particle  sizes,  each  assigned 
with  a  measured  statistical  frequency.  Usually,  the  effect  of  particle  size  distributions 
Is  small  if  the  distribution  is  narrow.  The  effect  of  a  large  fraction  of  larger  particles 
Is  to  act  in  the  reaction  zone  as  a  diluent.  An  experimental  and  analytical  study  of 
diluents  has  shown  that  such  systems  have  a  lower  reaction  temperature,  and  because 
of  corresponding  lower  diffusivities,  have  slower  reaction  rates.  Large  particles  have 
the  effect  of  increasing  ignition  times,  energies,  ami  delays.  Those  matters  will  be 
discussed  in  Section  IV. 


(3)  Effect  of  Phase  Changes 


The  computer  program  allows  for  a  phase  change  in  he  lower  melting  constituent. 
If  a  phase  change  occurs  in  a  system  of  stoichiometric  proportions,  no  effect  on  per¬ 
formance  is  observed.  If  the  phase  change  occurs  in  a  diluent  or  in  a  reagent  which  is 
present  in  excess,  a  phase  change  acts  against  the  reaction  by  lowering  the  reaction 
temperature. 

(4)  Effect  of  Changes  in  Composition 

As  illustrated  on  page  25,  changes  in  composition  are  handled  by  the  com¬ 
puter  program  through  the  selection  of  appropriate  diffusion  barrier  thick¬ 
nesses  which  arise  from  particle  sl2e  and  the  weight  fraction  of  the  individual 
constituents. 

(5)  Effect  of  Changes  in  Reaction  Mechanism  and  Properties 

Tliis  program  is  designed  to  describe  intermetallic  and  other  condensed  phase 
reactions  during  the  ignition  phase  and  assumes  that,  a  single  diffusion  rate  law'  can 
describe  the  process.  When  other  species  or  gas  phases  are  formed  as  a  result  of  the 
reaction,  variant  reaction  rates  must  be  expected  because  the  present  mathematical 
model  makes  simplifying  assumptions  as  to  the  constancy  of  transport  parameters, 
rates,  and  processes. 


SECTION  IV 


STUDY  OF  IGNITION  C HA RACT ERISTICS 


1 .  EXPERIMENT  A  L  STUDY 

Intermetallie  reactions  are  easily  and  reliably  ignited  by  brief  high-temperature 
events.  The  amount  of  ignition  energy  and  the  time  requirement  have  been  the  subject 
of  spei  ial  emphasis  in  this  program  because  of  the  great  importance  of  characterizing 
the  ignition  energies  in  many  types  of  pyrotechnic  processes.  Ignition  characteristics 
of  intermetallie  reactions  have  been  studied  experimentally  with  first  fires,  ignition 
wires,  spark  discharge,  and  thermal  flux.  The  common  observation  that  condensed 
phase  reactions  are  unreliable,  whether  thermitic  or  intermetallie,  has  been  shown 
to  be  wholly  caused  by  insufficient  attention  to  the  heat-transfer  characteristics. 
Systems  having  low  thermal  conductance,  such  as  those  containing  carbon  or  boron, 
have  been  shown  to  react  so  strongly  as  to  react  under  water,  once  initiated. 

In  Section  III  it  is  shown  that  excessive  thermal  conduction  causes  quenching 
of  the  reaction.  Similarly,  high  thermal  conduction  requires  high  ignition  energies. 
When  using  a  first  fire  or  hot  wire  to  ignite  an  intermetallie  reaction  mixture  such 
as  Ni-Al,  it  is  necessary  to  deposit  so  much  energy  into  the  mixture  that  the  tempera¬ 
ture  ol  the  entire  mass  is  raised  substantially  above  ambient  temperature.  Addition 
of  a  comjioneut  which  reduces  thermal  conduction,  such  as  boron  powder,  substantially 
reduces  the  amount  of  required  ignition  energy.  Low  thermal  conductivities  are 
therefore  one  of  the  requirements  for  promoting  ease  of  ignition.  Fur  this  reason, 
the  experimental  study  has  focused  on  the  boron  and  carbon  mixtures  of  titanium 
and  zirconium. 

Intermetallie  reaction  mixtures  of  titanium  or  zirconium  are  readily  initiated 
with  a  first  fire  which  may  either  be  pressed  as  a  layer  on  the  reaction  mixture  or  it 
may  be  poured  on  loosely.  Studies  were  made  using  A1/V205,  B/V2O5,  B/BaCrOq, 
Mg/BaCrOq,  and  composites  of  these.  The  amuunt  ol  lirst  fire  need  not  be  much  - 
large  masses  of  a  boron  and  carbon  mixture  with  titanium  and  zirconium  are  reliably 
ignited  with  a  small  fraction  of  a  gram  of  first  fire  which  is  triggered  by  a  hot  wire 
or  an  electric  match.  This  ignition  technique  was  not  developed  to  the  point  of 
obtaining  quantitative  data. 

Another  measure  of  the  ignition  energy'  was  obtained  by  stripping  the  glass 
casing  off  a  flashlight  bulb  and  immersing  the  tungsten  filament  in  the  reaction 
mixture.  Reactions  could  be  initiated  within  a  few  milliseconds  using  a  voltage 
which  depended  on  the  size  of  the  tungsten  filament  and  also  on  the  particle  size 
of  the  mixture.  For  2-p  boron  mixtures  with  titanium,  lb  to  20  V  would  suffice. 

Lower  voltages  would  cause  the  ignition  delay  to  be  increased  many  times  punting 
to  a  possible  change  in  ignition  mechanism.  High  voltages  of  about  100  V  would  tend 
to  shorten  the  ignition  delay.  The  current  and  voltage  input  to  the  filament  was 
recorded  on  an  oscillograph,  yielding  photographic  records  such  as  those  shown  in 


Figure  9.  The  ignition  time  (and,  hence,  the  ignition  energy)  was  assigned  to  the 
time  at  which  an  abrupt  change  in  current  and  voltage  was  recorded.  Such  changes 
were  ascribed  to  the  formation  of  a  thermally  conductive  reaction  mixture  or  the  burn¬ 
ing  through  of  the  filament.  This  method  does  not  necessarily  yield  the  minimum 
energy,  as  ignition  could  have  been  initiated  before  the  wire  burned  through  or  after¬ 
ward  due  to  the  continued  energy  deposition  into  the  reaction  mixture  from  the  power 
supply.  Ignition  characteristics  could  better  be  determined  through  controlled  timing 
of  energy  deposition,  or  by  other  techniques  such  as  laser  deposition  or  image  furnace 
technology.  Such  studies  were  beyond  the  scope  of  this  effort. 


o  20  40  f. o  ho  I  on 
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TIME  (msec) 


Ignition  Energy  Using  a  No.  328  Bulb 

18  V  x  0.80  A  =  14.4  W 

14.4  x  50  x  10-3  Bec  =  0.720  J 

*  0. 172  cal 

Figure  9.  Ignition  History  of  -325  Mesh  Titanium -Boron  Mixture 
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Am >t hut-  method  for  initiating  intermetallic  reactions  was  with  a  high  voltage  dis¬ 
charge  ,  A  feasibility  demonstration  was  made  by  discharging  a  2-pf  capacitance 
across  electrodes  which  were  inserted  into  a  line  titanium -boron  mixture.  It  was 
shown  that  the  liner  the  powder  mixtures,  the  lower  was  the  required  discharge 
voltage.  For  instance,  10-p  boron  was  initiated  with  the  capacitance  charged  to 
4000  V  while  a  60-p  titanium  mixture  was  not  initiated.  No  attempt  was  made  to 
monitor  the  actual  voltage  drop  across  the  electrodes  or  to  measure  the  current  of 
the  spark.  The  stored  energy  in  the  above  case  was  16  J  ,  with  the  capacitance 
discharging  in  0. 1  msec.  No  studies  were  performed  on  the  role  of  thermal  or 
electrical  conductivity  of  the  mixture  on  the  required  spark  energy. 

2.  ANALYTICAL  STUDY  OF  IGNITION  CHARACTERISTICS 

Ignition  energy  is  a  complex  function  oi  heat  capacity,  thermal  conductance,  and 
dif  fusivity.  It  can  be  expressed  as  energy  per  unit  mass  provided  that  the  mass  of  the 
(adiabatic ally)  reacting  material  is  specified,  or  as  energy  per  unit  area  provided  that 
the  ignition  time  is  also  specified.  Both  methods  are  illustrated  in  this  report.  The 
calculation  of  the  ignition  flux  and  ignition  time  by  a  mathematical  model  (p.  26)  is 
illustrated  in  Appendix  IV.  The  resulting  data  for  various  reaction  mixtures  are 
listed  in  Table  XV  and  shown  schematically  in  Figures  10  and  11.  No  exact  experi¬ 
mental  verification  of  these  calculations  is  available,  but  points  obtained  by  tungsten 
filament  method  (Figure  0)  and  by  ignition  in  a  filament  image  furnace  show  that  the 
computed  data  are  consistent  with  experimental  observation.  It  should  be  remembered 
that  the  ignition  fluxes  are  averaged  values,  explaining  why  the  most  reactive  niLxtures 
seemingly  require  the  highest  flux. 

In  examining  the  data  shown  in  Table  XV,  one  must  envision  variant  ignition 
times  and  fluxes  for  each  source  temperature.  This  means  that  those  data  are 
illustrative  only  and  they  should  be  repeated  for  each  ignition  method.  These  data 
yield  insight  into  the  magnitude  of  flux  required  and  so  facilitate  the  design  of  hardware 
or  of  verification  tests.  Both  ignition  time  and  fluonco  are  affected  by  particle  sine 
and  conductance,  explaining  the  results  obtained  in  the  high  voltage  discharge  study. 

Another  way  of  describing  ignition  energies  is  by  computing  the  mass  of  ignitor 
required  to  initiate  the  reaction.  As  before,  the  results  will  be  determined  by  conduct¬ 
ance,  d iflii sivity,  and  particle  size  as  well  as  by  the  temperature  (enthalpy)  of  the 
ignitor.  Tho  computer  program  for  performing  this  calculation  on  a  concentric 
cylindrical  geometry  is  shown  in  AppendLx  III,  and  it  is  illustrated  in  Figure  12. 

The  scheme  is  to  vary  the  size  of  the  igniter  until  the  resulting  enthalpy  is  sufficient 
to  initiate  a  self-sustained  reaction. 

Hither  method,  the  fluenee  calculation  or  the  ignition  energy  calculation,  may 
be  used  to  design  h  rdware  and  to  assess  the  safety  ot  pyrotechnic  mixtures  with 
respect  to  thermal  transients.  Other  geometries  may  be  similarly  programed,  and 
codes  can  be  written  which  calculate  ignition  energies  directly. 
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TABLE  XV.  COMPUTED  IGNITION  CHARACTtKlSTICS 


Contact  Time  With  Heat  Source:  1  sec 

Temperature  of  Heat  Source:  800  aK 

Thermal 

Conductivity 

Particle 

Size 

(M) 

Ignition 

Energy 

(cal/g) 

Ignition 

Time 

(msec) 

Ignition 

Fluence 

(cal/cm2) 

Ignition 

Flux 

(cal/cm2sec) 

D  =  1  exp  (- 

20000/KT) 

0.  00058 

3 

788 

5.  75 

0.90 

156 

0.0008 

10 

3710 

12 

1.79 

160 

30 

4130 

100 

5.  55 

55.5 

60 

4000 

464 

10.  85 

26.5 

0.0035 

3 

642 

7.65 

3.  06 

400 

s,» 

680 

24 

5.  75 

240 

0.  017 

10 

2280 

14.5 

6.  80 

470 

0.  054 

25 

5500 

740 

148 

200 

0.  16 

25 

5900 

370 

190 

515 

D  -  10  exp  ( 

-20000/RT) 

0.  0008 

1 

700 

0.21 

0.  01 

50 

0 

1640 

2 

0.  79 

400 

30 

1660 

19 

2.4 

126 

60 

1300 

74 

3.  72 

51 

0.0082 

38 

_ 

7770 

320 

48.  2 

150 

As  illustrated  in  Section  m,  computer  modeling  permits  the  assessment  of  the 
slative  importance  of  the  various  parameters  in  the  initiation  process.  Furthermore, 
it  may  serve  an  a  tool  for  studying  the  safety  of  pyrotechnic  mixtures  toward  ignition 
by  thermal  transients  or  spark  discharge. 

The  analytical  treatment  assumes  an  igniter  (or  first  fire)  of  known  density, 
heat  capacity,  thermal  conductivity,  and  which  is  brought  instantaneously  to  a  uniform 
ignition  temperature.  The  temperature  of  the  igniter  then  drops  uniformly  due  to 
heat  or  ..duction  into  the  interm  etallio.  The  ignition  energy  was  computed  for  an 
Igniter  of  L/D  -  „  which  defines  its  mass,  enthalpy,  and  surface  area. 
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REACTION  MIXTURE 


IGNITER 


a.  End  Burning  Ignitor 


b,  Internal  Igi.itor 


Figure  12.  Illustration  of  Mathematical  Approach  to  the  Determination 
of  Minimum  Ignition  Energy 


The  computed  ignition  energies  have  been  plotted  on  Figure  13  as  a  function 
of  thermal  ditfusivity  with  particle  size  as  a  parameter.  These  data  are  obtained  for 
diffusion  prefix.  D  ,  of  5  cnr/sec  and  an  activation  energy  of  20,000  cal/mole. 

Actually,  many  systems  have  lower  diffusivities  with  the  result  that  the  energies 
shown  in  Figure  13  would  be  expected  to  be  higher,  The  data  are  plotted  versus 
thermal  diliu.sivity,  a /p Cp  ,  in  order  to  compensate  for  variations  in  density. 

Figure  13  also  shows  experimental  data  for  ignition  energies  which  were 
determined  using  the  tungsten  filament  method  descv  oed  in  Section  IV.  The  agree¬ 
ment  between  experimental  and  analytical  data  must  be  considered  adequate  as  neither 
the  theoretical  diffusion  parameters  nor  the  actual  densities  of  the  test  mixtures  are 
known  precisely. 

The  main  observation  which  one  would  make  from  such  a  correlation  is  that  the 
required  ignition  energies  are  indeed  small  -  about  0.2  cal  or  1  J  for  a  typical 
titanium  mixture  with  carbon  or  boron.  This  corresponds  to  less  than  a  milligram 
of  the  typical  first  Ere,  or  a  diameter  of  an  igniter  cavity  of  less  than  a  millimeter. 

That  these  magnitudes  are  correct  is  borne  out  by  the  success  of  using  a  fine  tungsten 
wire  as  the  ignition  system . 

Such  findings  have  profound  implications  for  the  manner  in  which  conventional 
ignition  systems  are  used.  It  appears  that  greater  reliability  or  a  smaller  ignition 
system  would  be  possible  if  the  role  of  thermal  diffusivity  in  the  ignition  process  were 
better  appreciated.  For  instance,  conventional  thermite  systems  could  be  easily  ignited 
if  the  highly  conductive  aluminum  were  encapsulated  in  a  thermally  nonconduotive 
medium.  Recent  developments  have  shown  this  to  be  possible  (45).  Appendix  IV 
describes  how  the  computer  output  can  be  utilized  in  obtaining  the  ignition  energy 
flux  and  the  ignition  energy. 

3.  APPLICATION  OF  IGNITION  STUDIES  TO  INCENDIARIES 

The  preceding  section  described  experimental  and  theoretical  techniques  by 
which  the  effective  performance  of  intermetallic  reactions  as  igniters  could  be 
determined.  This  section  will  describe  briefly  several  demonstrations  of  these 
findings. 

a.  Ignition  and  Reaction  of  Thermites  and  Thennates 

Figure  14  shows  a  magnesium  cylinder  filled  with  Thermate  TH-3  and  fitted  with 
a  tungsten  filament  and  a  Ti-2B  igniter .  Voltage  signals  of  as  little  as  16  V  from  a 
dry  cell  will  ignite  the  system  within  milliseconds,  assuring  complete  combustion  of 
the  thermate.  This  type  o£  device  was  assembled  to  demonstrate  the  potential  for 
igniting  normally  sluggish  reaction  mixtures.  The  same  system  also  works  with 
thermite. 

Figure  15  shows  a  large  (250  g)  pressed  body  of  a  Ti-2B  mixture  into  which  a 
1/4-in.  cavity  was  drilled.  Into  this  cavity  was  inserted  a  tungsten  filament  Ti-2B  igniter 
assembly.  This  type  of  system  will  also  ignite  within  milliseconds  as  was  shown  by  a 
motion  picture  study.  The  product  cf  this  reaction  is  a  solid  clinker  with  a  fairly  uni¬ 
form  distribution  of  voids.  The  incendiary  action  can  be  augmented  by  dispersion  of 
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Figure  13.  Ignjtion  Energies  of  Intermetailic  Reaction  Mixtures 


Future  14.  Assembled  Incendiary  Demonstration  Unit 

1/2  x  2  in.  magnesium  cup  with  No.  32b  flash 
light  bulb.  Inserted  tungsten  filament  is  in 
contact  with  Ti-2I3  igniter  mixture,  the  remain¬ 
ing  cup  filled  with  therm  ate  TH-3. 


the  reacting  mass.  The  moat  effective  agent  for  breaking  up  large  chunks  of  inter- 
metallic  reactants  lias  been  the  addition  of  a  thermite.  An  example  is  a  mixture  of 
Ai/CuO  which  upon  r i>>e  generates  copper  vapor.  This  vapor  is  though),  to  be 
the  agent  responsible  for  the  explosive  action. 


The  high  coat  of  intermetallie  reactants  may  restrict  these  reaction  mixtures 
to  the  role  of  igniters  of  the  inexpensive  but  sluggish  thermites  and  the r mates.  As 
shown  in  the  preceding  section,  by  proper  control  of  the  thermal  transport  charac¬ 
teristics  of  '(Kith  intcnnetaiiic  and  thermite  mixtures,  greatly  improved  and  economi¬ 
cally  attractive  Incendiaries  can  be  designed. 

b.  Use  of  Intermetallie  Reactions  in  Tracer  and  Incendiary  Munitions 

Conventional  tracer  munitions  utilize  a  metallic  element  such  as  zirconium 
or  magnesium  with  an  oxidizer  such  as  a  chlorate,  perchlorate,  or  nitrate.  Problems 
relating  to  storability,  reliability,  and  operational  difficulties  have  prompted  a  search 
for  alternate  reagents. 

A  thermal  analysis  of  interrnetallie  reactions  showed  that  ignition  times  of  the 
order  of  dwell  time  in  the  rifle  (2.8  msec)  and  reaction  times  cm  the  order  of  the  flight 
time  (2  sec)  were  feasible.  Figure  10  shows  the  simulated  temperature  history  of  a 
tracer  bullet  which  was  initiated  in  5  msec  at  2800 “K.  Upon  emerging  from  the  barrel, 
the  temperature  of  the  exposed  surface  drops  while  the  heat  diffuses  into  the  bullet. 
After  approximately  40  msec,  the  temperature  has  risen  to  the  steady  state  value. 

This  shows  it  is  not  necessary  that  the  reaction  mixture  remain  exposed  to  the  heat 
source  until  the  reaction  temperature  has  been  reached  but  only  until  sufficient  beat  has 
been  added  to  overcome  the  diffusion  limited  reaction. 


Such  a  reaction  mixture  can  act  as  a  tracer  component  as  well  as  a  delay  fu/.e. 

Upon  impact,  a  visible  signature  may  result.  The  delay  fuze  may  ignite  a  pyrotechnic 
charge,  be  it  a  signal  flare,  a  smoke  grenade,  or  an  incendiary.  The  projectile  may 
be  a  bullet,  armor  piercing  projectile,  or  a  shaped  charge. 

Tests  were  performed  on  intermetallie  reaction  mixtures  which  were  uncompressed. 
These  will  react  with  the  formation  of  a  luminous  fountain.  The  use  of  such  reaction 
mixtures  as  filler  for  shaped  charges  has  been  suggested. 


The  feasibility  of  using  intennetallic  reaction  mixtures  as  igniters  under  transient 
conditions  was  demonstrated  by  the  insertion  of  u  pressed  interrnetallie  component  into 
a  smoke  cartridge.  Ignition  with  a  black  powder  charge  was  demonstrated  (40)  and 
resulted  in  reliable  generation  of  smoke  owing  to  the  retention  of  the  reaction  energy 
by  the  ignition  system  . 

c.  Formulation  of  Interrnetallie  Reaction  Mixtures 

Preliminary  tests  have  indicated  the  feasibility  of  formulating  and  shaping 
Interrnetallie  reaction  mixtures  in  a  continuous  manner.  Addition  of  glues  such  us 
gum  arable  and  Inorganic  binders  (45)  have  produced  firm  shapes  suitable  for  use  in 


igniters.  A  preliminary  test  has  shown  that  intermetallic  mixtures  can  be  extruded 
as  rods  or  tubes  when  these  are  mixed  as  a  paste  with  lubricants.  Fifty  percent 
by  weight  of  water  containing  5  percent  gum  arabic  or  a  similar  amount  of  volatile 
organic  solvent  containing  a  small  percentage  of  a  nonvolatile  binder  are  recommended 
for  this  purpose. 


APPLICATION  OF  INTERMETALUC  REACTION  TO  PYROTECHNIC 
AM)  ORDNANCE  TECHNOLOGY 


1.  CHARACTERISTICS  OF  INTKRMETALLIC  REACTIONS 

This  study  has  shown  intormetallic  reactions  to  have  many  characteristics  which 
are  desirable  for  pyrotechnic  and  ordnance  technology.  Several  such  characteristics 
have  been  measured  and  put  to  use  in  present  devices;  others  await  additional  develop¬ 
ment  effort. 

•  Shock  Resistance.  Mixtures  of  metal  powders  cannot  be  initiated  by  shock  as  by 
sudden  compression,  by  dropping  of  a  weight,  by  friction,  or  by  dynamic  shock, 
in  this  uav  they  differ  from  pyrotechnic  igniters  containing  chlorates  or  iron 
ordnance  initiators  containing  primary  explosives. 

•  Temperature  Res istance .  Mixtures  of  metal  powder  will  not  be  initiated  at  moderate 
temperatures  .2  This  means  that  pyrotechnic  components  containing  intermctallic 
reaction  systems  can  ho  stored  and  operated  at  temperatures  above  those  containing 
nitrates,  chlorates,  and  other  conventional  oxygenators.  Igniters  containing  inter- 
metallic  mixtures  can  operate  beyond  the  stability  limits  ol  ther mates. 

«  Corrosion  Resistance.  Py  rotechnic  and  incendiary  systems  containing  magnesium 
have  low  storability  because  of  the  effect  of  moisture.  Other  examples  are  the 
action  of  the  decomposition  products  of  single  and  double  base  propellants  or.  ordnance 
components.  Sensitive  materials  now  require  inhibitors  or  mechanical  barriers 
which  lower  reliability.  Intermctallic  reaction  mixtures  are  not  so  affected. 

•  Moisture  Resistance.  Intennetallie  reaction  mixtures  will  not  cake,  leach  out, 
diffuse  out,  or  reaet  under  the  influence  of  high  humidity,  moisture,  solvents,  or 
potting  compounds .  In  this  way,  they  ditfor  from  pyrotechnic  components  containing 
perchlorates  and  peroxides,  thereby  opening  a  new  field  lor  pyrotechnics 
development. 

•  Spark  Resistance.  While  intennetallie  reaction  mixtures  can  lie  initiated  by  a 
sizable  electrical  discharge,  they  are  not.  atfeeted  by  the  magnitude  ui  electrical 
energy  contained  in  stray  rf  signals,  or  static  discharge  by  operator  handling. 


“Reaction  mixtures  containing  the  alkali  metals,  mercury,  or  low  melting  alloys  of 
tin.  indium,  etc.  ,  are  exceptions.  As  shown  in  this  report,  high  thermal  conductivity 
mid  low  reaction  energies  make  these  systems  undesirable. 
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The  necessary  safety  requirements  inherent  in  handling  finely  divided  metal  powders 
must  be  met,  but  control  of  electrical  resistivity  of  the  reaction  mixture  can  be 
achieved  without  significant  change  in  operating  characteristics. 


•  Chemical  Resistance.  Intermetallic  mixtures  and  their  reactions  are  moderately 
insensitive  to  contaminants  and  variations  in  composition.  This  means  that  graded 
compositions  of  intermetallic  reaction  mixtures  with  other  pyrotechnic  components 
are  possible.  Low  quality  control  standards  with  resulting  low  cost  are  permissible 
for  many  applications. 

•  Condensed  Phase  Reactions.  Condensed  phase  reactions  are  not  affected  by  environ¬ 
mental  pressure  and  exposed  surface  area.  This  means  that  igniters  for  solid 
propellants  for  vacuum  ignition  can  be  made  and  that  underwater  application  will  not 
stop  intermetallic  reactions.  Lack  of  gaseous  products  causes  the  heat  of  reaction 
to  be  confined  to  a  small  volume  so  that  intermetallic  reactions  have  a  high  energy 
density. 

•  Compatibility  With  Various  Support  Media.  Intermetallic  reaction  mixtures  of  low 
thermal  conductivity  can  serve  as  delay  mixes  and  fuzes  without  being  affected  by 
heat  loss  through  the  container  or  mounting.  It  appears,  subject  to  further  testing, 
that  intermetallic  reactions  can  compete  with  existing  delay  mixtures,  but  operate 
more  compactly,  more  reliably,  and  in  a  more  adverse  environment  of  shock, 
corrosion,  or  temperature. 

•  Flexibility  in  Reaction  Characteristics.  Reaction  rates,  ignition  rates,  and  ignition 
times  can  be  varied  by  changes  in  mixture  characteristics.  This  means  that  a 
mixture  can  be  tailored  to  meet  specific  performance  requirements. 

2.  PRESENT  APPLICATIONS  OF  INTERMETALLIC  REACTIONS 

At  present,  intermetallic  reactions  are  chiefly  employed  for  the  ignition  of 
pyrotechnic  and  ordnance  devices. 

Pyrofuze  has  been  used  for  the  ignition  of  propellants  (4,  47);  also,  it  has  been 
studied  as  a  delay  time  fuze  (48), 

Nickel-aluminum  has  been  used  as  igniter-mixture  for  thermite  cutting  mixtures 

(49). 

Titanium-boron  inserts  have  been  tested  as  igniters  for  M84A1  HC  Smoke 
Cartridges  (46)  . 

3.  FUTURE  APPLICATION  OF  INTERMETALLIC  REACTIONS 

Bettor  understanding  of  the  means  of  controlling  reaction  characteristics  will 
bring  alxmt  a  greatly  expanded  range  of  applications. 


47 


i 


Ip.terineta.ilit;  reaction  mixtures  are  being  developed  as  active  components  in 
tracer  munitions  (50)  where  they  promise  by  absence  of  particulate  reaction  products 
to  prevent  jamming  of  the  machine  gun  mechanism.  Intermetallic-  reactions  when  user 
in  tracer  munitions  and  armor  piercing  munitions  promise  a  visible  signature  upon 
impact.  Intermetallic  reactions  have  been  mentioned  for  incendiary  projectiles  (45), 
as  low  ignition  energy  igniters  for  thermite  cutters,  as  infrared  heat  sources,  for  in 
situ  welding  of  space  structures  or  refractory  metals  (51),  as  shock  insensitive 
initiators  of  ordnance  boosters,  as  igniters  for  signal  markers,  as  active  components 
to  be  ejected  by  shaped  charges,  and  for  delay  fuzes  (48).  bimetallic  fibers  allow 
the  burning  rate  of  a  solid  propellant  to  be  increased  without  a  significant  decrease 
of  specific  impulse  (52). 

Use  of  intermetallic  reaction  mixtures  in  state-of-the-art  thermato  incendiaries 
will  assure  complete  and  rapid  ignition  with  high  reliability,  producing  an  incendiary 
which  is  economical  and  immune  to  countermeasures  such  as  physical  removal  from 
a  target  prior  to  ignition.  A  thermato  device  which  is  completely  ignited  will  react 
re|  ardless  of  the  condition  of  the  impact  area  which  may  lie  muddy,  sandy,  or  soft. 
Reduction  of  the  number  of  duds  will  improve  the  etfectiveness  of  present  incendiaries 
at  little  difference  in  cost. 

Intermetallic  reactions  are  completely  compatible  with  pyrotechnic  and  thermite 
reactions  which  tend  to  disperse  the  powdered  reaction  mixture  by  splattering,  Ignition 
times  of  thermite  and  thermato  systems  are  shortened,  and  their  rate  of  heat  output  is 
increased  so  that  improved  performance  in  visual  (signal)  and  incendiary  applications 
can  be  expected. 

Recent  improvements  in  Ihornnte-typo  reaction  mixtures  (45)  have  made  possible 
the  expansion  of  thermite  reactions  to  many  types  of  incendiary,  pyrotechnic,  and 
ballistic  munitions,  The  analy  tical  and  experimental  techniques  developed  in  this  pro¬ 
gram  are  directly  applicable  to  an  evaluation  of  other  pyrotechnic  mixtures  offering  a  now 
technique  tor  selecting  such  thermal  and  kinetic  properties  as  to  assure  sue;  essl'ul 
performance  of  the  devices  with  a  minimum  of  trial-and-error  testing. 
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CONCLUSIONS  AND  RECOMMENDATIONS 


Exothermic  intermetallic  reactions  f>  ve.  been  shown  to  be  an  altogether  new  class 
of  pyrotechnic  system.  Like  the  well  knov.  .1  e?  m.j;  systems,  they  react  without 
forming  a  vapor  phase.  The  reaction  r  Les  a:  e  determined  in  part  by  diffusivity  and 
particle  size.  Additional  parameters  >.hich  offer  a  possibility  of  control  are  thermal 
conductivity,  heat  of  reaction,  and  composition. 

Intermetallic  reactions  have  been  shown  to  be  amenable  to  successful  analysis 
of  the  heat  transfer  process  through  computer  programing.  This  has  been  proven  to 
be  a  powerful  tool  in  the  selection  of  components  which  will  produce  the  desired  reaction 
characteristics.  Moreover,  the  techniques  developed  promise  to  be  of  value  in  analyzing 
the  performance  of  other  condensed  phase  reactions  such  as  those  of  thermites. 

Of  particular  value  has  been  the  assessment  of  the  minimum  ignition  energy 
requirement  and  response  delay  characteristics  in  v  mis  of  particle  size,  thermal 
conductivity,  and  heat  of  reaction.  These  studies  ha\  e  pointed  the  way  to  a  wide  range 
of  applications,  ranging  from  utilization  of  into  •■rota’ lie  reactions  in  tracer  munitions, 
as  igniters  for  smoke  cartridges  and  signal  flares,  n  infrared  decoys,  to 
shock  and  temperature  resistant  boosters  for  ordnance  systems. 

Further  computer  studies  are  necessary  to  define  ignition  characteristics.  It  is 
recommended  that  computer-generated  findings  of  ignition  characteristics  be  corrobor¬ 
ated  by  experimental  study  of  intermetallic  as  veil  as  thermite-type  systems.  Such 
studies  should  be  performed  using  transient  energy  sources  such  as  spark  discharges, 
powder  flashes,  lasers,  and  shuttered  image  furnaces.  It  is  also  recommended  that 
further  pyrotechnic  and  ordnance  applications  of  intermetallic  reactions  be  tested, 
particularly  those  in  need  of  a  temperature-stable,  corrosion  and  moisture-resistant, 
and  shock-insensitive  ignition  system.  The  incendiary  characteristics  of  existing 
devices  can  be  improved  through  the  use  of  intermetallic  ignition  systems. 
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APPENDIX  I 

SOURCES  AND  SPECIFICATIONS  OF  REAGENTS 


Aluminum 

Alcoa  No.  120  Atomized  40%  -  325  mi  *h  15  q  Effective 

80%  -  -t0  i  >esh 

Alcan  MD  44  Atomized  21%  +  100  ms&c  30  <i  Effective  S?ze 

2C%  4  .200  mesh 
13%  +  325  mesh 
40%  -  325  mesh 

Nickel 

Sherritt  Gordon  Mines  Ltd . ,  Grade  F  5%  +  250  mesh  40  /u  Effective  Size 

40%  +  325  mesh 
60%  -  325  mesh 

International  Nickel  Co.  ,  Type  123  Carbonyl  Nickel  4-7  // 

Type  255  Carbonyl  Nickel  2.6  -  3.4// 

Titanium 

Atomergic  Chemical  Co.,  Approximately  10  </  size 
Consolidated  Astronautics  -30  +  100  mesh,  100  y,  Effective  Size 
Atlantic  Equipment  Engineers,  -50  *  100  mesh,  60  //  Effective  Size 

-200  +  325  mesh,  30  //  Effective  Size 
-325  mesh,  10  n 
10  q,  1  n  Effective  Size 

Zirconium 

Atomergic  Chemical  Co.,  Approximately  10  q  size 
Consolidated  Astronautics,  -30  +  100  mesh,  300  q  nominal  size 
Foote  Chemical  Co. ,  G  Grade,  2  q  size,  3  q  Effective  Size 

Carbon 

Graphite  No.  41B,  Consolidated  Astronautics,  0.5  q  Effective  Size 
Shaningan  Chemical  Co.,  Actylene  Black,  Agglomorates  1  q 

Fines  <  10-2  q 

Boron 

Research  Grade,  U.S.  Borax  and  Chemical  Co. ,  95%  Research,  Effective 
Size:  2  q 

Atlantic  Equipment  Engineers,  90—  92%  Grade,  Effective  Particle  Size:  10  q 
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APPENDIX  II 

PROGRAM  STATEMENTS  FOR  LINEAR  RATE  STUDIES 


Program  Description 

The  program  performs  the  following  operations  in  sequence: 

(1)  Calculation  of  the  temperature  in  each  length  increment  (slice)  of  a  rod- 
shaped  reaction  body  over  a  series  of  time  increments. 

(2)  Calculation  of  the  extent  of  reaction  in  each  slice  over  J  series  of  time 
increments. 

(3)  Use  of  a  measured  spectrum  of  particle  sizes  in  calculating  the  rate  of 
diffusion  and  the  rate  of  heat  evolution. 

(4)  Determination  of  the  temperature  rise  on  the  basis  of  data  on  phase 
transitions  for  the  low  melting  metal. 

Input  Parameters 

The  following  data  are  read  in: 

C  heat  capacity,  cal/deg-g 

RO  density 

KT  heat  conductivity,  ca,/cm-deg-sec 

DO  preexponential  factor  in  Arrhenius  expression 

E  activation  energy,  cai/mole 

Q  heat  of  reaction,  cai/g 

EPS  emissivity  (currently  set  at  0.9) 

R  radius  of  rod,  cm 

XL  length  of  rod,  cm 

XSOFO  nominal  thickness  of  low  melting  metal  in  mixture  particle,  cm 
XHARO  nominal  thickness  of  high  melting  metal  in  mixture  particle,  cm 
THKO  nominal  particle  size,  cm 

CONS  initial  concentration  of  low  melting  metal 

DX  thickness  of  each  slice  (distance  increment),  cm,  not  less  than  THKO 

FT  number  of  DDT  in  DT,  chosen  to  be  20  (see  discussion) 

TEMIN  temperature  of  ignition  source,  °K 
TIMIN  heating  duration,  sec 

TO  initial  temperature  of  sample  (usually  300°K) 

P  constant  (around  0. 10) ,  selected  to  adjust  magnitude  of  DT  (see  discussion) 

MXT  maximum  number  of  time  increments  desired 
TMLT  melting  point  of  low  melting  metal,  °K 

HFU  heat  of  fusion  of  low  melting  metal,  cal/cm^ 
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Discussion  of  Program 


The  distribution  of  particle  sizes  is  also  specified  in  the  input  as  follows.  The 
data  are  made  up  of  15  sets  of  3  numbers  for  15  size  increments.  The  first  number  is 
the  normalized  frequency  factor  for  the  ;  Lze  increment.  The  second  is  the  thickness 
of  the  low-melting  metal,  and  the  third  is  the  thickness  of  the  high-melting  metal 
in  each  particle.  In  this  manner,  variations  in  density  and  composition  of  the  mixture 
constituents  can  be  accounted  for. 

The  first  section  of  the  program  calculates  certain  quantities  needed  in  sub¬ 
sequent  computations,  sets  up  initial  conditions  for  temperature  and  for  the  unreacted 
material,  and  prints  out  the  input  data.  In  die  remaining  portion  of  the  program, 
a  complete  calculation  is  carried  out  for  each  time  increment,  DT  .  The  maximum 
number  of  these  increments  is  specified,  usually  1030. 

For  each  time  increment,  heat  conduction  is  first  calculated  and  a  tentative 
temperature  distribution  along  the  rod-shaped  sample  is  obtained.  Material  dif¬ 
fusion  is  calculated  independently.  This  is  the  most  time-ccnsuming  part  of  the 
program.  To  calculate  the  extent  of  diffusion  during  a  given  time  increment,  DT 
is  divided  into  a  number  of  smaller  increments,  DDT  .  That  number  is  determined 
by  the  program  for  each  individual  DT  so  that  DDT  is  small  enough  to  give  acceptable 
accuracy  in  the  diffusion  calculation,  yet  not  so  small  as  to  make  the  computer  time 
too  long.  Usually  the  number  of  DDT  in  each  DT  is  between  10  and  25. 

Since  the  rate  of  reaction  is  assumed  to  be  controlled  by  the  rate  of  diffusion, 
once  the  extent  of  diffusion  is  calculated,  it  can  be  used  directly  to  compute  the  heat 
evolved  and  the  fraction  of  material  that  is  still  unreacted. 

The  diffusion  calculation  is  carried  out  IS  times,  once  for  each  particle  size. 

The  heat  evolved  from  different  particle  sizes  is  summed,  using  spectral  coefficients 
as  weights.  Then  the  temperature  rise  due  to  reaction  is  computed. 

Thermal  radiation  is  also  calculated  independently.  It  is  assumed  that  during  each 
time  increment,  DT,  the  lateral  surface  of  the  slice  radiates  at  the  temperature  pre¬ 
vailing  at  the  beginning  of  DT.  The  heat  lost  by  radiation  is  distributed  uniformly 
throughout  the  volume  of  the  slice. 

The  tentative  temperature  is  now  corrected  for  heat  gain  from  reaction  and 
heat  loss  from  radiation. 

After  printing  out  the  results,  the  program  goes  to  the  next  time  increment  DT. 

The  magnitude  of  DX  is  important  for  the  success  of  the  run.  This  is  because 
the  time  increment  DT  is  calculated  from  it.  The  magnitude  of  DT  is  selected 
(through  the  equation:  DT  =  P  *  C-  *  RO  *  DX  **2/KT)  so  that  the  reaction  takes  place 
during  the  series  of  time  increments  specified.  If  DT  is  too  short,  the  reaction  will 
start  only  near  or  even  after  the  end  of  this  series.  The  velocity  of  the  reaction  front 


cannot  be  evaluated  properly  it  I)T  is  too  long  because  the  reason  will  be  completed 
in  all  slices  during  the  first  lew  time  increments,  In  this  case,  evaluation  of  the  speed 
of  the  reaction  front  might  not  be  possible  and  calculation  of  heat  transfer  and  material 
diffusion  might  be  too  inaccurate.  DX  is  chosen  to  be  approximately  equal  to  the  nominal 
particle  size. 


The  program  statements  follow. 


l  -• 

1* 

■V  - 

2* 

3* 

4  * 

£ 

9* 

^  # 

i  • 

4  * 

.<*  ■ 

Jljfc  • 

4* 

1C* 

ff. 

u* 

12* 

13* 

$ 

14* 

,ib 

19* 

** 

14* 

«■ '  * 

17* 

14* 

*/:<<- 

IS* 

,  -.ft 

2C» 

21* 

22* 

23* 

jW 

24* 

29* 

.jrf 

24* 

27* 

0 1  HP '  J  Of  J  I(100i3)i  0<100|3|i  F!l00»3* 
nlKPHSIOM  Spr<  15>«  XSfU5)>  XwR?15) 

Of  HP  1 1*!  ION  X?3F<15).  XH*R(1>).  TWKt  1*S )  •  S?CO< 13 > *  H0VSU»> 
plHPMSptl  T  HSOF  <  ).  00  ’  3 1 15  )  i  PftCX«  l00-:3<  i5  > 
p J MFlvIS !  Ofi  SOLF<300tl) 

REAL  KT 

300  FORMAT!!' F.10.3> 

4Q0  FcSu^J'c0’  .E7,2^X,6w  RO  •  J T  ..E7.2^X.6h 

*,.«  *.c7.2»?X.6h  £  s  »E7.2»?y**n  Q  -  *  P.7  •  4 » 2X  >  6hEPS  %?r 

402  KR’fcTUH  .lRX^H^*n!»S.K.^u»'lfiTH-«.HHICH  TEMP  <  5X ,  Hh1  GN  T,he.5X 

•  ,9Ho6qT  S]ZE<^Xie‘HnE.UT*x'  7X  i  6hdf.LT  A  7  > 

404  FCR'-uTtlH  ,EX#6HACTUAl<4Xi7lEiOi4i3X)//)  „  4X  7 

:  c“  rS* « w,  ^  ^  J  r:s; ;  :  Sht«So^,  D  ’ s 

wno.sx 

•  *  9  M7O0MTFMR  i  4X  i  4  HFT  INPUT  »6X.4H’1AXTI:-1> 

412  FCRMAT(1H  .rXi6HACTiULi4X.6(eiQ,4.3X) . tl4!  . 

420  FCR'1aT(  lH0.24X,4HSPC0ilc'x<4NX<iOr,l5X,4HXHAR»l5X.3H7HX> 

422  FCR;UT(  1H0»  ?OX,  4  <  f.10|4 , 10.X> ) 

430  FORMAT!  I  6 ) 

434  fCRmaTC  1M0*^x!iAHREI.T,NC  P'MNT  ■  .flU.l»2X.nHDEGRF.E  K'5X' 

»17hheaT  Of  FUSION  *  *F10j 1,2X>9hcaU7'cH**;>> 

RFAn<S,432>  (SPCtl>*  *sr<!>»  XHR(  ! )  t  1*1*15) 

KASr  •  0 

200  KASF  «  KASE+l 

1P'(KASF.  •(•<»*  1>  WR!TEt6i44U) 

440  FCRHAT(  1>|1»10hTiB2  30  MU) 

IF<KA$f  -F«,  2>  WRtTE<6,440 

442  FCR'4 4T(  1H1 » 10HT I  R2  6Q  MU 3 

RF  ao  < 5 * 300 )  Cl  Ro,  Kt  » t)0  i  c  O'  EPS 
RF  AO  d i 3  '0 )  R i  XL i  XSOFO.  XmAROi  ThKOi  CONS*  QX.  FT 
RFAo  (5.302)  TFM[N,  TJMIN.  TO*  P.  MXT 
RCA0<5.3  0)  t.nlt.  Hru 
c 

c  update  particle  size  data 

IF(kaS£  .EG,  1)  DEN  »  1333.3 
I  F  ( K  A  >F  •F.d.  ?)  DEN  =  64ft,6 
PC  i50  L  ■  1*  15 
SPCP'D  5  SPC<Di/24<’l 
XSOf  ( l  >  =  XSF<L)7’OEN 
XHAR<L>  2  XhR<  L )  7 0Efl 
150  CONTINUE 


nl”  tot.  NUiiapR  Op  Sfcti°n9  in  Ro0.  inclUdIng  *■  added 

VSECT  3  VOLUME 'of  SECTION 

aLaTS  =  LATERAL  SURFACE  ARE*  cf  section 

aenqs  *  r.NP  surface  area  of  sfction 

SGI-IF.PS  *  SIGMA.EPS1LON 

XL  *  length  of  rod 
r  *  radius  of  rod  „  „  , 

XSOF  ■  thicknfss  of  soft  mftal 

XHtO  *  THICKNESS  OF  HARD  metal 

CCNF  *  INITIAL  concentration  cf  soft  metal 
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51  * 

C 

TFMIN  *  10M(T1CS  TF.mPEAa 'j-<l 

62* 

C 

T  I M  J  N  2  1  ON  J  7 1  ON  7ME 

61® 

c 

MXT  s  MAXI M J’l  NUMBER  OF  TlMf  INCRE4FNtS 

61  * 

c 

SPen  «  COEFFICIENT  or  SIZE  DISTRIBUTION 

6  J» 

c 

ThkO  •  NOMINAL  PaSTICLE  S[Ze 

6J» 

c 

65® 

VSEcT  *  3.14i59»R®»2®OX 

64® 

4UATS  =  6 , 20 310  »R«OX 

67* 

AEWnS  =  3)14 159  *  R®  ®  2 

68  » 

SGMgPS  =  1,37E.12*FPS 

NL  *  IF  1  X<  XL^OX  > 

69® 

fa* 

IFUL'D*  "  FLOAT(NL)  iGTi  0,3!  Ml*NL+1 

71® 

NL  *  N|_*2 

72® 

MLM  =  NU'l 

73® 

00  136  1-1*15 

7®® 

THK(  1  >=XSOF( I >*XHAR! J < 

75® 

hovs<  t )*xhar<  n-'xsori :  > 

74* 

136 

CONTINUE 

77® 

PROVO  =  5.QE-O0 

78® 

qT  ■  P*C»RO®DX»*2/kT 

7  i  * 

ICTlf-l  =  1 F I X  <  Timim/dT> 

«o» 

IFdlMlN  .LT.  0T)  IqTIM  *  1 

81® 

HF'J  *  HFU®XSOf/<XSOF+VhAR! 

82® 

c 

01® 

c 

INITIALIZE  T( 1,1)  AMo  F( !iH  IN  THE  SECTIONS 

84* 

c 

F(  1 ,  1 )  iUn  F(N|_;1!  W 1 UL  NOT  gf  NEEDED 

83* 

1(1,1)  *  T  EM  I N 

84® 

DO  1  '0  I  =  2< NL 

87® 

T(  1,1)  •  TO 

85® 

F ( 1 , 1 )  «  1-0 

09® 

SCLF ( I , 1 >  «  1,0 

90® 

loo 

CONTINUE 

91* 

c 

92® 

c 

INITIALIZE  THSOF < I ■ 1 1 L )  AND  PROX(t,l,L>  IN  ThE  SECTIONS 

DO  1  '2  1  =  2, NLM 

93® 

94® 

DO  138  Ls 1  *  15 

95® 

TWSOF( I i 1,L>  *  XSOF<Ll 

94® 

PROX ( 1 ,1*L1  *  PROVO 

97® 

138 

CONTINUE 

90* 

102 

CONTINUE 

99® 

WRJTE  (6,400)  C*RO»KTtDO*E*Q«ePS 

100* 

WRITE  (6.402) 

101® 

WRjTEt 6,404)R,XL*TEMlN,TlMlN.THK0iDX*DT 

iol® 

WR l TE (6,410) 

103® 

WRITE* 6, 4 12)  XSOFO,XHARa«CONS,P,TO,FT*MXT 

104® 

WR J  T£( 6 , 4  34 )  TMLT,  HFU 

105* 

WR I TE( 6 , 4  20  5 

104® 

DO  140  L=1,15 

107® 

WRJTE* 6, 422)  SPCO< L ) i X$OF< L ! . XHAR< L> * THK* L> 

108® 

140 

CONTINUE 

109® 

0 

lie® 

c 

start  calculation,  thermal  conduct,  is  independent  of  temp 

111® 

c 

start  00  LOOP  FOR  TIME  JNOfcX  U 

112® 

DO  1  ‘6  JS1.MXT 

11.3* 

c 

m  and  n  Take  up  values  i#2  and  3  calculated  from  j 

114  * 

M  *  MOD*  U,3 )  *  1 

11?* 

N  «  M0D<U»1,3)*1 

114* 

c 

qO  uOOPs  FOR  SECT J0NS  Of  R°0 
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lt»* 

il«* 

11*- 

120-> 

121* 
122* 
123* 
124  * 
125* 
121* 
127* 
123* 
129> 
130* 
131* 
133* 
133* 
131* 
139* 
134* 
133* 
139* 
133* 
140* 
141* 
142* 
143* 
144* 
149* 
144* 
i  43* 
144* 
149* 
ISO* 
16<* 
152* 
153* 
154* 
159* 
153* 
152* 
158* 
159* 
160* 
161* 
162  * 
163* 
164< 
*69* 
166* 
16?* 
1 68  * 
169* 
170* 
171* 

172* 

172* 

174* 


C  DIFFUSION  COEFFICIENTS  CHANGE  WITH  J  *Ni)  t  BUT  NOT  WJTh  K 

DO  1  '8  !=2»NL* 

DO  148  l  ■  1 > 15 

IFITHSOFl -N£.  QsO>  UO  TO  U2 
THSOFt  1 . 1-1 . L >  »  0.0 
PROX<|«R*L>  -  THK<U 
GO  TO  140 

132  *  D0*EXP(-E/2«/TU »N) > 

KN  *  0 
DDT  =  DT/rT 
124  KN  *  KN+1 

DOT  =  DOT/KN 

TAM  =  DU,H)*0DT/PR0X<1»N|LJ 

!F(TaN  ,gT.  Oi 5*TnS0F i 1 <M|L )  .aNd«  DT/ddT  .LT.  25.)  G°  T0  134 
THSOF<I«M,L>  «  THSOFl  1 
PROX ! I »  N. l  i  =  PROXU,N,L> 

KFT  =  I  FI «<  DT/pDT ) 

c  material  diffusion  during  each  ot 

do  124  JL  «  2i  KFT 

Tam  *  D< I »M)*00T7PR0X< I »M|U 

THSOFl  I  iM’iL)  *  THSOF  l  I  »M,i.J-TaM 

IF(THSOF< ! *  M . U  >  -LT.  0 1 G  >  TwSOf(IiMiL)  *  0.0 

PROX  I  t  <  M.  l.)  =  PROXl  i*M>l'*TaM*I1.  *  nOVS<l_>> 

IP<PROX<  I » M.  U  3  iGT,  THKU>>  PRQX(I.M,L>  8  ThK<L> 

I F <  THSOF<  1 »  M .  t  >  .qT.  1 1 0EoO6*xS0e(  1.5  )  5O  TO  12* 

THSOFl I iMjL>  *  0.0 
PROX  I  | .  M.  L )  *  THK<L> 

GO  TO  1<8 
124  CONTINUE 

143  CONTINUE 
C 

C  CALCULATE.  FI  I  iM  )  AND  CHANGE  IN  F(!,M) 

3XS0F*0 ■ 7 
STHS0F=0. 6 
oc  142  L  «-  1-15 

STHSOP  »  SThSOf*SPcOID*ThSOF(  I  iM*L> 

SXSOF  »-  SXS0F*SPcO<L!*X3OfU.! 

1<2  CONTINUE 

Fl 1 1 M)  •STNSOF'SXSOF 
UELF  »  F<  f iM)-F I J 1 N ) 

c 

C  CALCULATE  heat  EVOLUTION  CCaL6Cm»«35 

HEV  ‘  -!>ELF*<RO*G  ♦  ll.O'SOLFl  I  »N>*HfU>  > 

C  CALCULat€  HEAr  rADIaT!ON  (CaL7cM*«3> 

IF<t  ,e0.  2'  GO  TO  U4 

I F i  1  .pG.  NLh>  GO  TO  ISO  • 

GO  TO  l13 

114  :F'J*!gTim)  118,118,120 

118  HHA  8  AUTS*SGMEPS*<T« I»N)*«4  -  TO**4) 

GO  TO  116 

120  HrA  8  I ALA”S  +  AENDS'»SGMEPS*IT( I • N 1  *  *  4  -  T0**4) 

116  HRA  *  HRA^OT/VSECT 
DELH  *  HEV  -  HRA 
C 

C  CALCULATE  HEAT  CONDUCTION  INDEPENDENTLY'  OF  RaDUTION  ANP  REACTION 

144  T  ( I » M  >  *  P»<T<l*l»N)*'T(J-l'N>)-(2.»P-l.’*T<l.N) 


1?**  ' 

176*  146 

l’r* 

17  9* 

179* 

1  SO®  164 

181* 

182®  162 

183* 

184* 

185* 

186* 

187* 

189* 

189*  186 

190* 

19 1 » 

192* 

193*  186 

194* 

193* 

194*  170 

197* 

194* 

199*  172 

2  00* 

201* 

202* 

203* 

204* 

2uS*  124 

206* 

207*  10B 

209*  C 
209*  c 
2lC* 

211* 

312* 

21,3*  C 
214*  c 
2l3» 

216* 

217* 

219* 

219* 

•error* 

220* 

221  * 

222*  C 
223*  126 

224* 

225* 

224* 

227* 

228* 

229*  128 

230*  106 


REASONS?  iT  JRL  FUR  F -C I  AT  S  0\  AND 

T(  1 1  f  | )  ■  T<  It*')  *  OE.L"</C/-*0 
TTAM  a  T ( l i V ) 

I P  t  r  <  r .  ^ >  *1T,  1  •  0E"O8  )  CO  TO  106 

I P  <  T  T AM*T ( J » H  )  )  162.164,166 

SCLF< ! . ^ >  *  SOLF<  J.M)  . 

CO  TO  108  'fQj. 

IF<  T<  I .  M )  .  LF.«  TMUT)  GO  TO  164 
HFM  =  (T<  !»'4)-TMLT>*c*R0 
HFMT  *  HFU*F< I iM>«SOIF< J i N) 

IF<HFN  i  G£ .  UFMT>  GO  TO  1*“ 

SCI.F c  !  •  M »  »  SO|_F<  I,N)»(  1,  O-hfm/hf^T) 

T<  I  ,M)  «  TMI.T 

GO  TO  108 

SCLF< !.M>  »  0,0 

HFM  =  HFlNftFMT 

T<  I .  M 1  *  TM|.T*’HFM7c/r0 

CO  TO  108 

IF<TT*M  -ST.  TMUT)  GO  TO  170 
SOUF<  !.M>  *  SOlf<  1 .  N ) 

GO  TO  108 

1F<  T  (  J.M)  .  |„T i  TH1T)  GO  TO  172 
SOL.F*  t » M  >  a  SOE.fr  6  J  *  N  7 


*0. 


CO  TO  108 

H«F.  =  (THLT*T<  I,M>5*c»«0 

H3UT  *  HFU*F<  J  .M>»<  l;0»SCLF-<  I  .N>  > 

IF<  HRE  • GE  >  HSUT)  GO  TO  174 
T < 1 1 M )  *  TMLT 

SOLF( t »  M )  *  SOlf<I.N)  *  ( 1<0-SOLF{ I'N) )*HRe7h9UT 
GO  TO  108 

T  ( 1 . M >  *  TMLT«(hRE-HSUT!/c/80 

SOUF< I.M’  *  1 ■ C 

CONTINUE 


BOUNDARY  CONDITIONS 
T ( 1 , M )  *  TEMI N 

l F ( J  ,GT.  IGTIM)  T( 1 , M 5  «  T(3,M> 
T(NL.M)  =  T<NU-2,M) 


for  printout  pattern 

I F  <  J  ,UE.  100  .AND.  MOD(JiIQ)  .EQ.  0)  GO  TO  126 
I F  <  J  . LE •  50  i ANQ,  M0d(Ji5>  R,  0)  qO  TO  126 
I F ( J  > UE •  10  lANQs  H0n<Jf  '  i  0)  GO  TO  126 
I F ( J  . E« •  1!  GO  TO  126 

1F(J  i |_E •  800  .AND.  MOO U.-  -EQ,  OJ  GO  T6  126 

the  non  arithmetic  statement  is  misspelled- 

IF<J  .QT.  800  .AND.  MOQ  ( J , 6Q  )  , eq ,  o)  GO  TO  126 
GO  TO  106 

HEADING  FOR  PRINTOUT  of  AhKAy5 
WRITE  (6.406) 

DO  124  1=2, NL* 

TIME  *  J*OT 

DIST  *  ( FLOAT ( I )  *  1,5)*DX 
II  1  l-l 

WR  I TE !  6 »  4G8  )  J,  TIME,  II«  01ST.  F  *  I » M » »  T<  I ,  M  >  *  D<I*M>»HRA 

CONTINUE 

CONTINUE 


231*  CO  TO  20" 

232*  END 
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APPENDIX  m 


PROGRAM  STATEMENTS  FOR  ANALYTICAL  IGNITION  ENERGY  STUDY 


Program  Description 

The  program  handles  a  cylindrical  geometry  where  the  ignition  source  is  a 
cylinder  concentric  with  a  cylindrical  reaction  mass. 

The  distance  increments  are  the  equal  thicknesses  of  concentric  shells  around 
the  igniter  while  the  radius  of  the  cylinder  is  the  counterpart  of  the  sample  length  of 
the  program  described  in  Appendix  EL 

The  igniter  may  consist  of  a  single  distance  increment.  If  so,  its  radius  and 
volume  can  be  varied  by  varying  DX.  Alternately,  the  igniter  may  comprise  a  specified 
number  of  cylindrical  shells  so  that  its  radius  can  be  varied  independently  of  DX. 

The  input  procedure  is  the  same  as  shown  in  Appendix  H.  The  form  of  output 
is  also  the  same. 

This  program  does  not  handle  the  problem  of  particle  size  distribution;  rather 
it  calculates  material  diffusion  on  the  basis  of  a  single  particle  size.  This  particle 
size  is  the  same  as  the  nominal  particle  size  given  in  the  input  for  the  other  program 
and  corresponds  to  the  most  probable  particle  size  in  the  measured  distribution. 

Phase  transition  is  not  taken  into  account  as  it  was  in  the  program  for  linear  rate 
studies  on  intermetallie  reactions. 

Thus,  the  list  o'  „v  data  for  this  program  does  not  include  the  distribution 
of  pai'ticle  sizes,  the  _  Ring  point  for  the  low  melting  metal,  or  the  heat  of  fusion. 

The  program  statements  follow. 
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mt 

% 


fi' 

0 

M 


'  # 


A 


1 


Nir-jcp-  Tcifj,3»,  0(333,3).  F  <  3^0 . 3  > 


1 


J* 

A. 

7* 

*• 

!• 

1C* 

n» 

12* 

I. 3* 
1«* 
15* 

14* 

17* 

II. 
If* 

2C« 
21  • 
£!• 
S3* 

24. 

s;« 

.34* 

S3* 

23* 

29* 

3C* 

31* 

32* 

33* 

34* 

35* 
34  * 
37* 
3«* 
3«* 
4C* 
41* 
42* 
43* 
44* 
4»* 
44* 
47* 
4  4* 
49* 
SC* 
51* 
>2* 
33* 

54. 

SJ* 

54* 

37* 

34* 

59* 


f  r  T  t  3C  '  .  3  ) .  1 4  T  (  3^0 ,  J  > 
fil-fiMTi  T-SOF  t  3 JO  .  3 > «  <  30?. 3) 

rf-'F'.p  i <vj  v ( i? , 

REAL  * T .  LA"1'1*.  <Tt 
3o;  rcH-'*Ti  ’r'liM 
3^2  fcr^ati  «r  lv.  3,  •  i j  1 

403  FC-i-itf  lul.f  .|  ’c  •  ,F7,2,?X,6b  PQ  >  ,-t.2.2*.6h  KT  « • E7 • ? • 2X • 6m 
*pC  « . r 7. ?, 2* • 4m  c  =  iE7<S«rXi6i.  (j  *  * F 7, 2» 2* • frH£PS  * 

432  fCH"4T(1.4  .  1  **  .  I^RanI  'S«  7X  •  *>m.E  iGtH»6Xi  9m!Gn  TfcNF' .  5X .  *h  1  G‘j  LNGT.5X 
• . 9wP4RT  S  ]  2c  <  5X,6nni  tT*X  i  7  4  i  6ur.FL.3  *T  ) 

4ll4  FC««/.T(Ui  .£X.6m*CTUaL»4**  7<E10.4,3X)//> 

406  pCH'*iT(///l^.?X,  /^T  j*.cfX,4x.7HoL  TflR.4X.7nX  |  NCEX  ,  4X  ,  7^  qJST. 
•3*  .  7()U'JPr  4C*F.3  •  JX »  7nMt  Tf'^tSv  .*ho<  I  .  "4  > .  6X. 3hDUT  .  8X.  XHTHSf’P, 

•  6X , <nP»pX/ ) 

4nS  FCR^TCIH  .  •'V.  ]4 . 5X.£7.4, 4X,  J4, 1x.8(2X.g9. 4)  ) 

4 1 C  PC»sATMH  ,  1M  ,  *HSOPt*V  ,  7X  .  4HH*'«L4T  .6X,  (\HSQF  C0NS.5X.HhR  RATIO, 3X 
•,8H|GN  rOMU.4X.RnPT  f.p!jT,6x.6H '*»XTM> 

412  PCH->*T(lil  .4X,6HACTlJ*L»4X.6tEl0i4i3X),  114) 

416  pCH-'AT( //10X  ,  11^  LTt  MP  w  .E9.416H«LTFriH,i4x,llw0LTIf4E  *  »fc9,4, 

•  6HRLTl!-'r.l«X.llHLLDlfT  1  ,  RP  .  4 . 4h«i.U  JR7  > 

4 18  FCR‘'AT(//l04>33i.r,Tf.T*/C;TA'J  •  i)2Tu Ta/UZF.7 *2  -  .  2X . E« ,  4 , 12MUF /&T 

«AL  -  .F.®-4.  13h'(  TET*1«TtT‘J4  )) 

430  F(*R'14T<  J  6  ) 

pA  T A  V<  1  )  .  V<  21  .  V<  3) .  V(  4  ),wt  5)/'  t  In0»  ’  .  MX.  ' • '  Jm’.'ii  «  »  *  .  '  13. 
»'  / 

DATA  v<4’»V<81,v<  10>.V<12)/'4x.  '  ,  *  C  1  >4A  1  •  '  »  'tlHXl.*.  Mim(j))'/ 
REAR  15.3po)  C.  RO.  K  T » [j  C .  E  «  0.  EPS 
RE *0  15.300)  R.  XL.  XSOF.  XH4B,  CONS,  DX.  FT 
RF4"'  <5,3"2)  TEMIN,  X 1  ,  XT  1  ,  P,  MXT 
C  DE Fill! NO  THE  CONSTANTS 

C  ALWAYS  SET  CPS  «  O.J 

c  R  *  TOT,  RAOfJS,  1'lf.LUDINr.  ICNlTfR  •  Xl*Xu 

c  ML  .  TOT.  N.i-lagB  ap  SHELLS  IN  CvLlNr)F.R.  INCLUDING  1  IN  IC^Te* 

C  XL  •  CTLTN[=!CaL  Thic<UESS  (F  S7MP|_E  NOT  INCLUDING  1 GN I T  f.R 

C  X5  0F  »  Twj  rXNp.SS  OF  SOFT  f'f  TAL 
C  XWaR  ■  Thj  C41ESS  OF  HARO  vtTAl 

C  CONI  »  initial  concentration  CF  SOFT  HCTaL 
C  TEMfj  i  Initial  T ENPfR aTUh t  Of  IgNITeR 
c  xi  ■  ratios  of  igniter 

C  KT  I  »  HFAT  CONDUCTIVITY  OF  IGMTER 

c  '-'XT  =  ilAXpun  NUMBER  OF  T  jNp  INCREMENTS 

TC  •  3 JU  •  n 
NL  •  I F  I  X  (  X  c/  OX  ) 

)F(vl'd*  •  fLOaT(NL)  ,GTf  0,5)  ''L«NL+1 

ML  *  NL  *  1 

MLM  c  NL* 1 

THK  *  XR0c*vhah 

HCVR  ■  XiijR/XSOF 

XPcRS  1  Twk/110, 

pP(jxo  1  StCE-OW 

DT  a  p» C*RO*nx.»2/i<T 

oxi  *  sortIcit.kti/IP'C.RO1 

AM  «  Xl/DXl 
F4CT1  *  1.  ♦  l./AM 
F*CT2  *  2,  *  1./4M 
*  IN  *  (XI  -  0X1 ) **2 
ACU  =  0X1»( ? . *XI "0X1 ) 

TEM°F  a  2,/e 


l 
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60  TIMFF  =  ?,*Q*D07<C*t*XS0F**2) 

61  •  DiSTF  *  S'5KT(2.«fi»oQ*RO/(KT*E*XSOF*i«2)> 

62*  GAMMA  *  Z.^C/E 

63*  L*MPA  *  0,0 

64»  c  INITIALISE  T  ( I ,  i  5  *N0  Ft!*1'  IN  The  SECTtONS 

65*  c  Fti.U  aMo  F<NU»1)  will  NO1  BE  ''EEdEd 

6 «»  no  100  !=1*MU 

bi»  Ttj.n  .  to 

68*  IF<  I  .£3.  1>  T< 1 . 1 J  •  TtHtN 

69*  IOC  CONTINUE 

7C*  Ft  1 • 1 )  •  0.0 

71*  pO  1~2  1=2. AIM 

73*  F< 1*1>  *  1.0 

73*  THSOF<  i»l»  «  *SOF 

74*  PROV< 1 >1 >  *  P«OXO 

75*  102  CONTINUE 

76*  WRITE  <6. 400)  c.RO<KTioO»E'0.pPS 

77*  WRjTg  <6.40?' 

7<t*  UR  j  Tp  (6.404  )  r?»  Xe<  TIMIN'*  Xp  ThK*  0X*  dt 

79*  UR1TE<6,410) 

8C*  UR  i  Tf; «  6  *  4 12  >  XSOp*  XhAR*  C°NS.  P.  KTj*  f-*T.  MXT 

hi*  WR 1 Tp t  6. 416  )  TEMpf,  TJHeF*  OiSTF 

«2*  WRITE* 6.418)  GA"Ha.  lahda  „  .  .  „ 

kj*  p  start  calculation,  thermal  cOnouct*  i9  indepence.nt  of  temp. 
84*  c  start  do  loop  for  time  index  j 

85*  nO  106  J=1«MXT 

oa*  c  h  and  h  take  up  Values  i»2  and  3  calculated  from  j 

87*  M  ■  M0ptJ.3>  *  1 

88*  N  ■  MOD* J-l*3)*l 

89 »  c  pc  loops  for  sections  cf  hod 

9C*  c  DIFFUSION  COEFFICIENTS  CHAnGE  WITH  J  aND  I  BUT  NOT  WITH  K 

9l«  DO  108  1=1. NLM 

92*  I  Ft  1  .EO.  1»  GO  TO  114 

92*  IF<  ThSOF<  I»  N)  .A'E*  0,0)  GO  TO  132 

94*  TMSOF< 1. M)  ■  O.Q 

93*  PROX<  I.M)  *  THK 

96*  GO  TO  110 

97*  112  0 *  I . M 1  *  D0*EXP<-E/2./Ttt.N>' 

94*  KN  »  0 

99*  DOT( I , M>  »  OT/FT 

IOC*  134  KN  *  KN* 1 

101*  DOT< I  *  M>  •  ODT< I.MI/KN 

102*  TAM  =  o{ 1 .M)«OOT' ! .M)/PROKt I -N) 

103*  1 F  <  TAM  .GT.  0,5*THSOF ( S .N) '  GC  To  134 

104*  THSOF< I . M>  ■  THSOFfl.N) 

lf)3*  PROX<  I .  M  '  »  PROX<  1  ,N) 

106*  JFKI.M)  •  IFlX<DT/DDTt  I«M' > 

107*  KFT  =  1  FT<  I  »fl) 

108*  c  MATERIAL  DIFFUSION  DUPING  Each  dt 
1Q4*  DO  124  L=2.<FT 

lie*  Tam  *  pt 1 »M)*ODTt I (M)KPROXt |  .F ) 

HI*  IF<TaM  ,GT.  THSOFM.M')  TA^*TwSOFf  1»M' 

U3*  TPSOFtl.m  »  TH50F<  IiM)-TAM 

H3*  PROX<  I .  M '  ■  PROXd.M)+TAM«>l.  *  H»VS  > 

114*  IFtTHSOFt I. Ml  ,5T«  l.OE'Oi^XSCF'  GO  TO  124 

115*  ThSOF<  I . II '  »  0.0 

11«*  PROX<  t,M>  ■  THK 

117*  GO  TO  11 ■ 


* 


m 


ff  n ' 

.*A 


ff 


t*-  V 


its* 

119* 

12C* 
121* 
122* 
123* 
12<* 
129  » 
126* 
127* 
126* 
129* 
13C* 
131  * 
132* 
133* 
13<* 
135* 
134* 
13  7* 
13S* 

1 3  9 » 
14C* 
1  4  1  * 
142* 

1 4  J  * 

1  4  4  o 
1  ■»  5  » 
J46* 
147* 
146* 
149* 
ICC* 
1M» 
lr>?  » 
153* 
lri4 » 
1*3!* 
134  » 
137  * 
134  • 
1S9* 
1*C* 
161* 
162* 
1  3  « 
164* 
169* 
14.4* 
167* 
lftfl* 
169* 
1  7  C* 
171* 
172* 
177* 
174. 
175* 


X24  CONTINUE 

C  CALCULAT(:  F  <  I  *  M  )  A  NO  CH4NGk  IK  F  <  1  »  M  ) 

Un  F*IiM>  a  TmSOE<  1  .M/XSGF 
OELF  e  pCj.MI-f* ltM> 

1 F !  nFLF  -P.T.  0.0)  DELE  ■  0*0 
C  CALCULATE  HEAt  EVOLUTION  <  CAL>CH**3> 

HE  V  =  -RO*G*[>eLF 

C  CALCULATE  HELr  RaOIaTION  T C aL/Cm**3 > 

DFLui  *  HEV 

c  calculate  hf4t  cONivictjon  !hdf.;,fndem-lv  of  raoIatIUN  Ar;-  Reaction 

IF(]  1 0£  •  2)  r-0  To  u A 

114  T  {  I  . :  l  >  ->  P»n  (J  ,H)+FACTl«Tlt+l.nn>tp*FAcT2-l.l*T(l.M) 

TAM  =  ( Y< I ; *AOU+T< I . M)*AiM'/Xl » *2 
T  i  i  ,  M  '  *  '  a1-' 

gC  to  ion 

116  T (  J  s il >  «  p*<  T(  ; -i,M  >+rACTl«TC  1  +  1,-J)  )-(P»FaCt2-1. '»T5  1  *N> 

C  RE  AOJUE  7  Tt--'F'FHATU«F  FOR  BADlAT|OH  AMD  EEACTI0l'l 
T  C  I . ; !  A  a  T \  t . M  >  +  OEL^/C'00 

105  CCNT 1  HUE 

C  ;9C')m[iaRV  CONDITIONS 

T  t  Ni.  i  3 )  =  TdlL-2.'n 

[  F  (  J  .  LC  •  100  .aNO.  hOp<J,10)  .  EC3  •  0)  CO  TO  126 
jr(.j  ,(_c.  3]  .  a  M[j  i  NOptj,S>  .£<3.  0!  CO  TO  126 
IE(J  *  L.F  *  1C  vAN;.«  U0p(j*2>  *e«i  0>  G0  TO  126 
1F(J  .pO-  1)  GO  TO  12P 
C  FOR  PRINTOUT  PATT£HN 

I F  ( J  ,LF.  sno  .AMO.  MPOlJt^O)  .r,<5.  O)  GO  TO  126 
IFtj  .GT-  50C  .ANO.  MOp ( J 1 1 00 7  iEO.  0)  GO  TO  126 
GC  TC  106 

C  he  Am  MG  FOR  PR]'., ’TOUT  Op  AR«AVS 
126  WRITE  (3.406) 

DC  12*  I =1 • VLN 
T]iif  *  J*OT 

D 1  ST  a  ( FLOaT ( ] )  -  1,9:#DX 

WR]Te<6.40»)  J.  TIME*  J)  UlST,  F  f I  *  rl ) •  T(]iM!»  pfjitfJ,  DDT  < I » M  > • 
•TMSOK  I.(').°P0X(  I  *  H  > 

12B  CONTINUE 

[)C  14?  ]--9,nLM 

n  *  t-i 

I E ( 1 1  ■ NF ,  1  .AND.  M00<1 1,10)  .Np,  0)  go  TO  1*2 
LSOr  a  !Fi:<;th10F<  I,Ml/XPf.HS) 

IF(THSOF(  ],M)/VFt.RS  »  FL0AT(L50F)  .GT.  0-5)  LSOF  *  LSUF*1 

IPtLSOF  ■ LE <  0)  L?0F°1 

LPKOX  =  ]FIX(RR0X<  I,M)/).PEHGi 

I  F  t  FPROXt  1  ,l-*)7XptRS  -  FLOATClPRoX)  .C.T.  0  *  13  7  L^OX  *  LPROX+1 

1  F  !  LPPOX  « LE  •  0)  LPr<UK,l 

LHA°  s  HO  -  LSOf  -  Lp«OX 

IF<LHaR  'LE'  0>  LHA«"1 

ENcC[)£{  pld  NT7 , 430  )  lS°F 

V<7  )  a  PRjmT7 

ENCODE ( PRIMT9, 430 )  LP*OX 

V<  9)  *  PRji'IT? 

EArppE<PRllillli430)  LHAR 
V<  :.l)  *  PR  JM1 
WRITE* 6.V>  II 
1<2  CONTINUE 

106  CONTINUE 
E.ND 


vjfw-‘ 
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APPENDIX  IV 


METHOD  FOR  COMPUTING  IGNITION  CHARACTERISTICS 


The  method  for  computing  the  ignition  characteristics  from  a  computer  run  is 
illustrated  for  a  Ti-C  mixture.  The  computer  output  has  been  plotted  on  Figure  8. 
Table  IV-1  lists  the  temperatures  at  the  center  of  the  first  slice  of  a  rod-shaped 
reaction  mixture  for  various  time  increments  until  the  temperature  of  the  first  slice 
equals  the  igniter  temperature  of  800°K.  At  this  time,  the  heat  flux  into  the  first 
slice  vanishes,  and  the  total  increase  of  sensible  heat  in  the  slice  is  the  ignition 
energy,  and  the  average  of  the  temperature  gradients  multiplied  by  the  thermal 
conductance  becomes  the  average  ignition  flux. 


TABLE  IV-1  METHOD  OF  IGNITION  FLUX  CALCULATION 
Input  Parameters 


Thermal  Conductivity:  0.017  cal/cm-sec-°K 
Heat  Capacity  at  800 °K:  0.20  cal/g-°K 
Half  Thickness  of  Slice,  1/2  Ax  =  0.004  cm 
Temperature  of  Source:  800°K 
Contact  Time  with  Source:  Continuous 
Unit  Time  Increment,  At  =  0.1742  msec 


Slice  Temp. 

T  (°K) 

Temp.  Diff. 
AT  (°K) 

Time 

(msec) 

Temp.  Gradient 
AT/1/2  AX 
(°K/cmXl04) 

Ignition  Flux 
(cal/cm2  sec) 

Change  of  Temp. 
With  Time 
(AT/AtxlO6) 

300 

500 

0 

12.4 

2120 

2.86 

341.7 

458.3 

0.1742 

11  5 

1952 

2.63 

376.4 

423.6 

0.3484 

10.6 

1800 

2.43 

430.6 

369.4 

0.6968 

9.24 

1570 

2.12 

470.2 

329.1 

1.0452 

8.24 

1400 

1.89 

527.5 

272.5 

1.742 

6.80 

1156 

1.56 

607.6 

192.4 

3.484 

4.80 

816 

1.10 

654.0 

146.0 

5.226 

3.64 

620 

0.838 

686.3 

113.7 

6.968 

2.82 

480 

0.650 

732.0 

68.0 

10.452 

1.70 

290 

0.390 

800 

0 

14.5 

G 

0 

0 

Results 


Ignition  Time  * .  14.5  msec  2 

Average  Ignition  Flux  over  14.  5  msec .  470  cal/cm  sec 

Ignition  Energy  ** . . . . . . .  2260  cal/g 


♦Ignition  time  is  less  than  ignition  delay,  as  ignition  delay  is  the  time  required  to 
approach  reaction  temperature. 

♦♦Energy:  fc 
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